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INTRODUCTION 


The contract covered by this report involved a number of separate 

but closely related activities related to the field of microwave radiom- 

etry. These activities included the following specific tasks: 

(1) The development of several control loop and dynamic noise model 
computer programs for simulating microwave radiometer measurements 
has been essentially completed. This effort was initiated under 
an earlier task order contract, and the results of this develop- 
ment have been published (refs. 1 and 2). 

(2) Computer modeling of an existing stepped frequency radiometer 
has been done in an effort to determine its optimum operational 
characteristics . 

(3) The classical second-order analog control loop has been investi- 
gated to determine its relative and possible optimum performance in 
reducing the error of the estimate in a microwave radiometer based 
on noise bandwidth and settling time criteria. Such results were 
not found anywhere in the literature. 

(4) Several designs of digital signal-processing units for a microwave 
radiometer have been investigated extensively. One particular 
design has been identified as the most promising, and its behavior 
has been simulated extensively by means of the digital simulation 
programs mentioned in (1) . 

(5) Efforts have been initiated for developing the hardware and software 
required in the implementation of the digital signal-processing unit. 

(6) Some of the general characteristics and peculiarities of digital 
processing of noiselike microwave radiometer signals have been 
investigated. 

(7) Technical support has been provided in the form of computer data 
reduction of results obtained in actual ice missions. 

During the period of performance, the following specific persons were 

involved at different times (all on a part-time basis) : William D. 

Stanley (principal investigator, ODU) , Roland W. Lawrence (graduate research 



assistant, ODU) , Sally E. Kerpelman (undergraduate research assistant, 
ODU) , John M. Jeffords (graduate research assistant, ODU), and William 
H. Thornton (faculty consultant, ODU). 

Dr. Stanley was directly involved in all tasks except (7), and 
Lawrence participated in all but (3) and (7) . Jeffords performed most 
of the work in (3), and the appendix of this report was authored by him. 
Prof. Thornton provided consultative assistance in tasks (4) and (5) . 
Kerpelman provided virtually all the assistance required in task (7) . 

This report provides a documentation of the efforts of tasks (2) , 
(3) , (4) , (5) , and (6) . An attempt has been made to divide the report 
into logical units commensurate with the specific concepts investigated. 
However, due to a significant amount of overlap in some of these areas, 
the report does not necessarily proceed in an exact chronological format 

Use of trade names or names of manufacturers in this report does 
not constitute an official endorsement of such products or manufacturers 
either expressed or implied, by the National Aeronautics and Space 
Administration. 
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CONTROL ANALYSIS OF FEEDBACK RADIOMETER 


As a prelude to the modeling, analysis, and design of a Dicke 
noise-injection feedback radiometer, it is desirable to delineate 
the particular manner in which the control mechanism functions. A 
simplified block diagram of the important parameters of the closed- 
loop system is shown in figure 1. The signals appearing on the figure 
represent only the "dc-like" control levels expressed in terms of 
the effective temperature values at different points in the system. 

The fluctuation components of the signals do not appear in this model. 

The model shown in figure 1 will hereafter be designated as a "control 
model." The more elaborate statistical model will be referred to as 
a "dynamic noise model." 

This simplified control model only shows one Dicke switch since 
it assumes that all frequency variations between the first Dicke switch 
and the correlation switch are sufficiently high that they can be 
neglected for control- loop analysis. The gain between the two Dicke 
switches has been momentarily lumped with the transfer function G(s). 

The frequency variation associated with G(s) is primarily that of the 
integrator and loop filter which constitute the estimation circuit 
portion of the system. 

Analysis of the loop will now be made. The input temperature T^ 
undergoes a very slight attenuation due to the transmission coefficient 
1 - of the directional coupler (typically 0.99). The injected noise 
is weighted by the coupling coefficient of the directional coupler 

(typically 0.01) and added to the antenna temperature. The effect of 
the Dicke switch is to introduce an additional 1/2 factor so that T^ 
on figure 1 can be expressed as 

Te = I (Ti - Tg) (1) 

where Tj is the effective output temperature of the directional coupler 
and Tg is the ambient temperature of the constant temperature enclosure. 
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Effective output temperature Ti in turn can be expressed as 



( 2 ) 


where Tj^ is the injected excess noise appearing at the attenuator output. 
In turn, Tj^ can be expressed as 


T' = 
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where is the attenuator constant, T„„ is the excess temperature 

AI 1 Ea 

of the noise diode, t is the pulse width of injected noise pulses, 
and Ki is a constant relating output voltage to pulse frequency 

F^. The output voltage transform can be expressed as the 

product of the transfer function G(s) and the error temperature trans- 
form T (s) . For convenience, the s-domain notation will be omitted in 
most of the subsequent calculations. Thus, 


V = -G T„ (4) 

o E 

Equations (1), (2), (3), and (4) may be solved simultaneously to 
yield a composite expression for the output voltage as a function of the 
input and reference temperatures. After some manipulation, the result 
is obtained as 

Vo = (i - 'Sdc ) [r ;% - 6/ - 2 - ] (^B - \) 

where, for convenience, 3 has been defined as 


= *'dc *^att '^ex 


( 6 ) 


The output pulse frequency F^ is related to the input and reference tem- 
peratures by 


F = 
o 



G Kj/2 1 
1 + G B/2_. 



( 7 ) 
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One other possible quantity of interest is the "error" temperature 
degined in terms of a hypothetical point following the addition of 
noise-injection temperature, but after the 1/2 factor of the Dicke 
switch. This function is 

T = 1 ~ ^Dc) (^A ' ’^b ) 

E 2 1 + G e/2 


E* 

the 


( 8 ) 


As far as the control mechanism is concerned, the significant frequency 
dependency of the loop is contained in the G(s) transfer function. 

A typical arrangement involves a near-perfect integrator and a loop 
filter such as 


G(s) 


K 


o 

sCl + sT) 


( 9 ) 


The presence of this particular function results in a type 1 control 
loop. The resulting steady-state error due to a step function change in 
temperature at the input is zero on a deterministic basis. On the other 
hand, the presence of the noise fluctuations and a finite input bandwidth 
prevent the statistical error from ever reaching a zero value. 



MODEL OF ANALOG RADIOMETER SYSTEM 


The particular radiometer used as the reference for virtually all 
the analysis and design work involved in this study was the 4.5 to 7.2 
GHz stepped frequency system developed by Harrington et al. (ref. 3). A 
simplified block diagram of the system is shown in figure 2. This 
particular radiometer is a Dicke square-wave correlated type employing 
closed- loop noise-injection feedback to reduce the effects of gain 
fluctuations. Many of the salient operational features of the system 
have been adequately described (ref. 3) and will not be discussed here. 
The emphasis in this development will be directed toward the establish- 
ment and utilization of certain mathematical models for describing 
various aspects of system operation. 

While much of the emphasis in this overall study was devoted to 
the development of potential digital signal-processing schemes for 
radiometers, it was apparent at the very outset that the operation and 
limitations of the analog system had to be fully investigated first. 
Consequently, considerable effort was expended in obtaining various 
transfer functions and operational parameters for components within 
the analog system. Each of several sections will be discussed 
individually. 

RF Subsystem 

The RF subsystem will be defined as that portion of the radiometer 
from the output of the directional coupler to the input of the square- 
law detector. For the dc control model, all frequencies are sufficiently 
high that the net effect of this entire portion of the system may be 
represented by a single constant gain factor. On the other hand, for 
the dynamic random model, it is necessary to utilize an appropriate 
transfer function to represent the effect of the wideband filter pre- 
ceding the square-law detector. 

An approximate gain constant for the RF subsystem is determined by 
the procedure that follows. During the time that early testing on the 
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radiometer was performed, R. F. Harrington of NASA/LaRC accumulated a 
significant amount of data relative to signal and noise levels used in 
the evaluation. Based on a specified noise figure of 5 dB, Harrington 
estimated that the effective receiver temperature was about T = 627 K 
at some nominal frequency in the operating range. The Dicke reference 
temperature was established as Tg = 308 K. With the reference tempera- 
ture applied continuously to the RF subsystem, the signal level at the 
input to the square- law detector was measured at approximately -33 dBm. 

The input power density producing this output is estimated to be 
k(TB + T^) = 1.38 X 10’23 x (308 + 627) = 1.29 x 10'20 W/Hz or 
-168.89 dBm/Hz. The corresponding gain x bandwidth factor is then 
estimated to be -33 dBm - (-168.89 dBm/Hz) = 135.89 dB-Hz. This 
corresponds to an absolute power gain x bandwidth of 3.884 x 10^^ Hz. 

As it turns out, an early calculation based only on T employed an 
input power density of -170.63 dBm/Hz, which resulted in a gain 
bandwidth factor of 137.63 dB-Hz or about 5.792 x 10^^ Hz. This value was 
rounded to 5.8 x 10^^ Hz and used in many of the system simulations. 

This difference of about two dB is within the uncertainty range of the 
actual gain factors of the various components in the forward part of the 
loop as well as the frequency dependency range, and is believed to be 
insignificant in view of the general uncertainty concerning overall 
forward gain system parameters. This point of view is further reinforced 
by the fact that the gain involved is in the forward loop and has a 
relatively small effect on the closed- loop steady-state accuracy of the 
measurement within reasonable limits. 

Diode Detector 

The diode detector employed is an HP 5082-2350 Schottky diode which 
is biased in a square- law region of operation over the entire operating 
range of the radiometer. Within this region, the output dc voltage is 
then a linear function of the input power. Measurements by Harrington 
indicated a 200-yV output when the input was about -33 dBm (or 0.5 yW) . 

The corresponding gain constant is thus 200 x 10*^/0.5 x 10“^ = 400 V/W. 
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Video Amplifier 

The video amplifier is that section of the radiometer between the 
square- law detector and the output loop filter. This amplifier is 
basically a low-pass amplifier, but it need not pass dc since the mean 
temperature estimate being measured appears as a modulated component 
on the Dicke frequency. However, for square-wave correlation, it is 
necessary that a relatively large number of components of the square- 
wave be transmitted by the video amplifier with minimum attenuation 
and phase shift. Furthermore, noise injection with short pulses 
mandates that the amplifier have a large bandwidth compared with the 
reciprocal of the injection pulse width. (This point is pursued further 
later in this report.) 

The actual circuits used in the video amplifier will now be con- 
sidered. The first two stages, designated as ICI and IC2, are identical, 
and each has the form shown in (a) of figure 3. The transfer function 
is determined to be 

SR 2 C 2 ^1 + sCi (Rq + Ri)] 

' Cl + sRiCi) (1 + SR 2 C 2 ) 

With the element values shown on the figure, the function is 


G (s) = 
a 


66 X 10"3s(l + 2.3892s) 


(1 + 66 X 10"3s) Cl + 79.2 X 10"^s) 


-3. 


(11a) 


66 X 10“ 3s [1 + s/C2ir x 0.0666)] 

[1 + s/(27r x 2.41)] [1 + s/(27r x 2.0095)] 


(11b) 


The corresponding Bode plot for ICI and IC2 is shown in (b) of figure 3, 
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The circuit diagram for stage IC3 is illustrated in (a) of figure 4. 
The resulting transfer function can be determined as 

sR2Ci(1 + R4/R3) 

G (s) = (12) 

1 + s(RiCi + R 2 C 2 + ^2^1^ S^RiR2CiC2 


With the element values shown on the figure and the variable resistor 
adjusted for maximum gain, the function is 


G^(s) 


10.5s 

(1 + 0.53047s) (1 + 28.2767 x 10"^s) 


(13a) 


^ 10.5s 

[1 + s/(2tt X 0.3)] [1 + s/(2tt X 5.63 X 10^)] (13b) 

The corresponding Bode plot is shown in (b) of figure 4. 

Combining the effects of two sections of ICl and one section of IC2, 
the complete Bode plot of the original video amplifier is shown in 
figure 5. Observe that a zero of transmission occurs at dc and that a 
low-frequency rolloff of 18 dB/octave appears. Furthermore, the response 
is relatively flat from above a few hertz to about 5 or 6 KHz. However, 
it was later discovered that the bandwidth should be broadened considerably 
to accommodate the relatively narrow noise- inject ion pulses. 

DC Estimating Circuit 

Following the correlation switch, the signal is processed by an 
estimating circuit. The estimating network originally consisted of an 
analog integrator and a loop filter, which was a lag- lead network. The 
lag- lead network was later redesigned as a straight lag network, and this 
process will be discussed in a later section. 



The circuit diagram of the integrator is shown in figure 6, and the 
circuit diagram of the original loop filter is shown in figure 7. The 
integrator has a transfer function given by 


G,Cs) = - 


1 

RCs 


-17,857 

s 


(14) 


for the element values shown. 

The transfer function of the lag- lead network is 


G Cs) = 


Ri + Rf 


1 + 


sRqRi C 
Rn + Ri 


1 + sRqC 


(15) 


For the particular element values shown on the figure, the transfer 
function is determined to be 


11 (1 + 2.545 X 10 3s) 

(1 + 28 X 10-3s) 

11 [l + s/(2tt X 62. 525)] 
[l + s/(2ir X 5.684)] 


(16a) 


(16b) 


The Bode plot of the lag- lead network is shown in (b) of figure 7. A 
Bode plot representing the combined effects of the integrator and lag- 
lead network is shown in figure 8. 

Noise- Inj ection Scheme 

The noise- inj ection process in the closed- loop system is achieved by 
the use of a voltage-to-frequency (V/F) converter, which produces constant 
width pulses whose frequency is directly proportional to the loop output 
voltage. The pulses are applied to a PIN diode driver circuit which is 
used as a gate for an avalanche noise diode. When the pulses are on, the 
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noise diode output is directly applied through an attenuator to the 
auxiliary coupling port of a directional coupler. Conversively, when 
the pulses are off, no injection noise is applied, and only the 
ambient temperature appears as an input from the feedback path. 

Under stable ideal closed- loop conditions, the injection noise plus the 
input antenna noise is just sufficient to balance the reference noise, 
and the effects of fluctuations in the forward path are eliminated. 

If the parameters of the feedback path are known accurately, the 
antenna temperature may be measured by determining the number of pulses 
used per unit time to establish the balance. 

Data on the noise diode indicates that the excess noise temperature 
is about 31 dB above 290 K, corresponding to an excess temperature of 
approximately 365,000 K. The V/F converter is a Datel System VFV-lOK, 
which has a sensitivity of about 863 pulses/s/V. Each pulse has a width 
of 70 \iS so that the constant relating duty cycle to output loop 
voltage is about 863 x 70 x 10”^ = 0.0604. The injected noise encounters 
a 6-dB fixed attenuator corresponding to a power gain constant of 0.25. 

The injected noise temperature adds to the ambient temperature of 308 K 
at the attenuator output. The resulting signal temperature then adds to 
the antenna temperature through the directional coupler coupling coefficient 
of -20 dB, corresponding to a power coupling constant of 0.01. 

The manner in which the various constants contribute to the feedback 
noise is illustrated in figure 9. The possible input antenna temperature 
range is assumed to be from 0 to 300 K. The corresponding ranges of 
different temperatures are illustrated on the figure. The duty cycle 
ranges from about 0.334 to 0.00868, and the corresponding pulse rate ranges 
from 4774 pulses/s to about 124 pulses/s. Finally, the output voltage 
range is from about 5.532 V to 0.143 V. The injected noise has a 
sensitivity of about 25.55 K/pulse at the diode itself or about 0.0639 
K/pulse at the point of injection. This corresponds to a sensitivity of 
about 22,050 K/V at the diode or about 55.124 K/V at the injection point. 
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The manner in which the pulses add to the unsymmetrical open- loop 
signal in order to balance out the net temperature is illustrated in 
figure 10. The waveform in (a) depicts the initial situation in which 
the reference temperature plus the receiver temperature is greater 
than the input brightness temperature plus the receiver temperature. 
However, the effect of added noise-injection pulses brings about a 
net balance of energy as illustrated in (b) before filtering and 
smoothing occurs. The levels in the tables are the typical values 
appearing at the input to the estimation circuit; and of course, the 
results are smoothed by that portion of the system. The fluctuations 
about the mean levels are not shown in these illustrations; i.e., 
these are the "dc-like" levels. 

Other components in the system whose effects have been considered 
as constant for the purpose at hand include a tunnel diode amplifier 
manufactured by AERTECH having a noise figure of 6 dB and a gain of 
26 dB ± 1 dB, a mixer and preamplifier stage produced by RHG having a 
noise figure of 7.5 dB ± 0.5 dB and a gain of 47 dB, and bandpass filters 
produced by CIR-Q-TEL having an insertion loss of less than 1 dB in the 
passband. 

Combining all the constants discussed in this section, a control 
loop model for the stepped frequency radiometer is shown in figure 11 
based on maximum gain for the variable gain stage. After several steps of 
manipulation, the closed-loop transfer function for this system is obtained 
from standard control loop analysis as 




VoCsD 

TgCs) - T^(s) 


20.214(1 + 2.545 x IQ-^s) 
s^ + 38.5496s + 1114.1 


(17) 


The poles of the transfer function are readily determined to be 
s = -19.27 ± j 27.25. This means that the damping factor is a = 19.27 and 
the damped radian frequency is = 27.25 rad/s or f^ = a)^/27r = 4.34 Hz. 
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The undamped natural radian frequency is = 33.38 rad/s and = 
to^/2ir = 5.31 Hz. The damping ratio is C = ~ 0.577. 

It can be shown that the equivalent one-sided noise bandwidth 
for this transfer function is given by 


rl + 

" [ 8 ^ 


( 18 ) 


where T 2 is the numerator time constant in the lag-lead network (2.545 ms 
in this case). The value = 7.28 Hz was calculated for the parameters 
just considered. 

The sensitivity AT of the closed loop feedback radiometer can be 
determined from the relationship 


AT = 2(Tg + T^) 



(19) 


where B^^^ is the equivalent output noise bandwidth of the loop and 
Bgj is the equivalent statistical bandwidth of the input wideband filter. 
Using values obtained from the system, the values of AT for the four 
possible input bandwidths are summarized as follows: 


Bgj. 20 MHz 100 MHz 500 MHz 2 GHz 

AT 1.60 K 0.71 K 0.32 K 0.16 K 

These values represent sensitivities based only on the noise-reduction 
filtering and do not include variations in coupling factor, stability of 
the noise diode temperature, or Dicke ripple, which will be discussed 
later. 
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REDESIGN OF THE ANALOG LOOP 


Although a primary objective of this study was the development of 
techniques for digital processing of radiometer signals, certain 
immediate operational objectives suggested that some design modifications 
of the existing stepped frequency analog radiometer would be worthwhile. 

The first area of concern was that of the loop filter following the 
integrator in the estimation circuit. Prior to this study, the circuit 
used for this purpose was the combination of circuits shown in figures 
6 and 7. This circuit is a lag- lead network whose transfer function was 
given in equation (16) . 

The lag- lead network, in contrast to the lag network, has a larger 
equivalent noise bandwidth for a given damping ratio and thus does not, 
in general, smooth the fluctuations as effectively. (Although, for the 
values initially used in this system, the difference was not appreciable) . 

A separate study was made by a graduate student, John M. Jeffords, to 
determine various tradeoffs and possible optimum strategies for designing 
a second-order loop in a feedback radiometer. While much standard 
information concerning the usual control design of a second-order loop 
is widely available, peculiar aspects of the radiometer design necessitated 
a separate study for this purpose. The results of Jeffords’ study are 
given in the Appendix. 

Jeffords used the equivalent noise bandwidth times the settling time 
as a criterion for optimization. This quantity will hereafter be referred 
to as the B^,t product. Jeffords found that the absolute minimum B,.t 
product did result from the simple lag network, and this occurred for 
5 «0.95. When a numerator lead time constant (denoted as T 2 ) is used, 
the constant increases rapidly. Simultaneously, the value of c 

which results in minimum B.,t also increases. 

N s 

An immediate reaction to this result might be to set T 2 = 0 corresponding 
to the lag network. However, other criteria must be considered before 
this conclusion is completely obvious, such as the sensitivities of c 
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and the steady-state error to changes in loop parameters. After all, 
the loop parameters, particularly forward loop parameters, are subject 
to fluctuations, so the stability of the error to changes in parameters 
should be considered. 


Jeffords found that the sensitivity of the damping factor as a 
function of gain variations is a minimum in the vicinity of C = 1 for 
03 ^X 2 =1. It was further found that the sensitivity of the error to 
changes in loop gain is smaller over a wide range of time for ^^^2 
products of 1 to 1.5 or so. The result of these last two findings could 
suggest the use of a lead network for certain applications in which a 
wide variation of gain fluctuations is expected. The results of Jeffords' 
study are developed in depth in the appendix. 

As it turns out, the variation in gain parameters for a typical 
radiometer is not expected to be excessive. Furthermore, if some addi- 
tional time above the minimum settling time is provided, the increase 
in error is not expected to be significant even though the sensitivity 
is high. Typical calculations suggest that this conclusion is correct. 

After considering various aspects of the loop properties and the 
results of the study, a decision was made to select a value of c in the 
vicinity of unity (i.e., a critically damped loop) using a lag network 
only. The critically damped loop has certain design advantages, and 
it is sufficiently close in behavior to the optimum value that most of 
the properties are essentially identical. The product for the 

critically damped loop is approximately 1.5. However, an additional 
time should be allowed for the loop to fully settle because of the 
higher sensitivity of the error to parameter changes in the case of the 
lag network. 


The redesign of the loop based on a critically damped response was 
then formulated as follows; If the lag- lead network is replaced with a 
simple lag network, its transfer function 


G/s) 


K 




1 + sT 1 


( 20 ) 
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The closed-loop transfer function is then determined to be of the 
form 


v„Cs) 

TgCs) - T^Cs) 


Kl 


s *^LOOP 

+ ■ 4 * 

Ti Ti 


( 21 ) 


where Kj is the forward gain constant and ^'"lqqp gain 

constant . 

Following a review of the settling time specifications, it was 
decided to design around a damping constant a = 50, thus producing a 
transfer function of the form 


Vq(s) 


Kl 

(s + 50)2 


( 22 ) 


Comparing equations (21) and (22), it can be deduced that Tj = 0.01 s and 
= 8.82. The circuit form used to realize the appropriate lag network 
is shown in figure 12. This circuit replaces the original lag-lead 
network. 

In order to properly study the effects on the actual radiometer of 
the lag-lead and lag networks, the test circuit of figure 13 was connected 
in the system, and the parameters adjusted to visually determine what 
appeared to be an optimum response. The results compared favorably with 
those obtained from the mathematical analysis. 

One other design consideration for the loop will now be discussed. 

The injected noise pulses appearing at the radiometer input are relatively 
narrow, i.e., they are 70-ys wide. In the Dicke correlation process, it 
is necessary that all noise injected into the first half-cycle actually 
appear in the first half-cycle at the output switch. Any energy that 
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"spills over" due to bandwidth limitations will result in some error in 
the measurement process. Because o£ the short pulse widths, it was 
discovered that the bandwidth of the video amplifier was insufficient 
to process the noise-injection pulses without some distortions. 
Consequently, a degradation in the calibration factor as a function of 
the noise level was observed. 

As a first step, the capacitor was removed. When the problem was 
still observed, it was noted that the particular operational amplifier 
used had an insufficient bandwidth. Consequently, it was replaced 
with an operational amplifier having a much higher bandwidth and slew 
rate, and the problem was thus resolved. 


DICKE RIPPLE ANALYSIS 


Along with the noise fluctuations present in the output of a 
Dicke radiometer, the ripple corresponding to the Dicke switching 
frequency and its harmonics must be considered. This ripple is present 
due to the modulation and demodulation processes of the switched signal. 
While these components may be kept to a sufficiently low level, they 
are often ignored in the design process, and this can lead to troublesome 
operation. Judging from an extensive survey of the literature, it 
appears that many investigators are not particularly concerned with 
their presence since virtually no discussion of their effects has been 
found by this investigator. 

For the purpose of the ripple analysis, consider the block diagram 
shown in (a) of figure 14. It is assumed that steady-state conditions 
exist, in which case the square-ware is perfectly balanced with a 
sufficient amount of noise injection at the input. The constant 
represents all gain up to the point of the Dicke correlator, and this 
gain is assumed to be frequency independent as far as the frequency range 
of the Dicke switch and its harmonics are concerned. Note that this 
constant does not include the 0.5 Dicke constant shown in figure 11, 
since that constant arises in the switching process and does not appear 
explicitly in the figure. 

The function G(s) represents the net transfer function of the 
estimation circuit, which includes the integrator and the lag-lead or 
lag network. All of the significant frequency dependency in the loop is 
contained in this block. 

The constant K 2 represents the feedback constant. As in the case 
of Ki, no significant frequency- dependent effects are assumed in this 
block. 

The signal v^(t) appearing at correlator output is shown in (b) of 
figure 14 under balanced conditions. Fluctuations about the positive 
and negative mean levels are not shown. Since the function shown is a 
simple square-wave, the harmonic components may be readily determined. 
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The fundamental sine or cosine component of a symmetrical square-wave 
such as this is 4/n times the peak amplitude. This conq)onent is then 
weighted by the transfer function evaluated at s = where = 

2Trf^ and f^ is the Dicke switching frequency. Letting Vi represent the 
peak output component at the Dicke frequency, this value is 

Vi = I Ki(Tg * Tj^)|G2(Jm_j) I (1 - (23) 

While the actual magnitude of the fundamental component is certainly 
important, a more meaningful measure is the relative level of the fundamental 
compared with the dc level. The dc level in the steady state is 

(Tr - Tj) 

= K2 » - Sc’ 

Let aj = represent the ratio of the fundamental ripple level to the 

dc component. This quantity is determined from equations (23) and (24) 
as 


«1 


Vi 

v~ 

o 


4 CS 


n (T, 


- T^) 

- S’ 


Kl|G2(jw^)|K2 


(25) 


The quantity Kj |G 2 (ju)^) IK 2 in equation (25) represents the net 
loop gain for that particular frequency without the 0.5 Dicke switching 
factor. In the dc control loop analysis, it was convenient and proper 
to include the 0.5 factor in the definition of the loop gain. In order 
to make this analysis consistent with the control model, the 0.5 factor 
is absorbed in the loop gain definition, and an additional factor of 
2 is put in front to compensate. Thus, the ripple factor aj can be 
expressed as 


otl 


8 

n (Tg - T^) 


Gt (jWj) 


(26) 
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where [Gj^(jo)^)| is the magnitude o£ the loop gain transfer function 
evaluated at the Dicke frequency with the 0.5 factor included. 

A symmetrical square-wave of the form shown in (b) of figure 14 has 

only odd harmonics. The ratio of the nth harmonic magnitude to 

the dc component V will be denoted as a . This function can be 
^ o n 

determined as 


a 


n 


8 

nn 


(T. + T ) 

CTb”- T^) I 


(27) 


For a typical loop gain function, the fundamental component is 

usually much larger than any of the other components. Consequently, the 
assumption that the total harmonic level is approximately the same as 
is usually valid. This assumption will be made in the analysis that 
follows. In the existing radiometer, the Dicke frequency f^ is 
approximately 120 Hz. 

For the original loop design, the ripple level of the fundamental 
can be determined as 


ai = — 


^ V 


51.195 (1 + 2.545 x 10~3jg)) 
jo)(l + 28 X 10"3j(^) 


X 100% 


(28) 


For the modified design, the corresponding function is 


ai 
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jw(l + O.Oljw) 


X 100% 


(29) 


As it turns out, the ripple levels for these two functions are almost the 
same, so only one will be shown graphically. The output peak ripple 
level expressed as a percentage of the output signal reading is shown 
in figure 15. This does not include the effect of post-loop filtering, 
which further reduces the ripple level. 
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From this curve, it is readily observed that the ripple level 
increases markedly with increasing input temperature and could be most 
troublesome if not properly eliminated. Computer simulation of the 
two systems did indeed verify that the ripple level was dominated 
by the fundamental component, and (b) equal to the value predicted by 
the mathematical model. 

One way in which the ripple level could be reduced significantly 
in the analog system is by means of additional loop filtering. In order 
not to disturb the control mechanism of the loop, which is dominated by 
two very low frequency poles, it is necessary that the loop filter 
have a cutoff or filtering range well above a few hertz (the range of 
the dominant poles) , but at the same time it must have a high attenua- 
tion at f , = 120 Hz. 
d 

Several possible analog filters were investigated for this purpose 
and simulated with the available computer programs. Some three-pole 
Butterworth filter designs with cutoff frequencies below 100 Hz were 
used, and the simulations verified that they would be quite effective 
in reducing this ripple. 

One particular filter of interest for this purpose was an analog 
notch filter whose notch was placed at the Dicke frequency. The 
particular filter used had the transfer function Gj^(s) given by 


- 


(1 + 1.7591 X IQ-^s^) 


1 + 2.5678 X lO'^s + 5.0555 x 10”^s^ + 4.375 x 10 ^s^ 


or 


^ ^ 402.124(s 2 + 568.49 x lO^) 

(j (sj - ■ 

^ s^ + 1.1557 X lO^s^ + 587.01 x lO^s + 228.60 x 10^ 


(30) 


(31) 


The amplitude response of this filter is shown in figure 16. This 
filter was found to to quite effective since it completely eliminates the 
fundamental component. However, there is a small effect due to the 
higher harmonics which are not attenuated by the notch filter. 
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GENERAL DIGITAL CONSIDERATIONS 


In this section several general considerations concerning digital 
processing of radiometer data will be investigated. Most of these 
considerations are general in that they apply to systems other than the 
radiometer, but it was necessary to review their significance in the 
context of the radiometer analysis. Some of these factors arise because 
of conversion between the analog signal and the digital signal and the 
reverse process. 

First, consider the analog system shown in figure 17 with an input 
x(t), an output y(t), and a transfer function G(s], where G(s) = 
Y(s)/X(s). Assume now that the input x(t) is chosen to be a wideband 
noise source; i.e. x(t) = n^(t). Assume that the ’’power" spectral 
density of is uniform and given by G^^(f) = V^/Hz (or W/Hz) 

on a one-sided basis, and assume that = 0. The variance of the 

output = E can be expressed as 


■ ■ 00 
y^(t) = E[y2(t)] =f |G(j(jo)|^n^ df 
o 


(32) 


This result can be expressed in countour integral form as 


y^(t) = 


■ml 


G(s)G(-s) ds 


(33) 


where C is a contour encompassing the jto-axis and an infinite semi- 
circle in the left-hand half-plane. It is further assumed that all 
poles of G(s) are located in the left-hand half-plane and that the 
high-frequency rolloff rate of lc(jw) | is no less than 6 dB/octave. 

Now assume that a sampled-data version of x(t) is formed, and let 
x(n) denote this quantity. Assume that all the samples of x(n) are 
statistically independent. (More will be said about this assumption in 
a later section.) Assume also that a discrete transfer function H(z) 
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corresponding in some sense to G(s) is formed, and let yCn) represent 

the discrete output. From the properties of z-transforms, Y(z) = H(z)X(z), 

However, before involving some of the properties of the discrete-time 

signal, consider the process of transforming the s-plane contour integral 

of equation (34) into the z-plane in a direct correspondence sense, 

assuming of course that H(z) corresponds closely to G(s). Since 
sT 

z = e is the relationship between the z and s variables, dz = 
sT 

Te ds = Tzds and ds = dz/Tz. Substituting this expression in equation 
(33), and replacing G(s) by H(z) and G(-s) by H(z“^), the following 
form of the contour integral results: 


y^(n) = 




(34) 


where y^(n) is the output variance of the sampled signal and Ci is 
the corresponding contour in the z-plane, 

A result from discrete transfer function theory is the relationship 


y^(n) = ^ 


/ 

L c 


H(z)H(z-‘) 




dz I (n) 


(35) 


for the case where all samples of x(n) are statistically independent . 

From the preceding results, it can be deduced that the quantity 
represents a power spectral density in V^/Hz for a continuous analog 
signal . However, n^/2T = (n^/2)f^ represents a total variance in for 
a sampled-data version of an analog signal, and this corresponds directly 
to x^(n). Thus, the relationships correspond to each other with power 
density (or voltage-squared density) being the input for the analog case 
and total variance being the input for the discrete case. In both cases, 
the output variance is the total variance in volts -squared. Looking 
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at it from a different point of view, the total input in is 
divided over the range f^ = 1/T for the discrete-time signal case. 

In systems involving A/D and D/A conversions, it is necessary to convert 
back and forth through these relationships to envision the manner in 
which the signal forms change. 

The point of this discussion is that the discrete-time and continuous- 
time forms produce identical results when viewed from the proper 
perspective. When continuous-time systems are considered, it is more 
convenient to work with power voltage-squared) spectral density 
and to consider the power divided over the frequency range of the 
sampling frequency (or half the sampling frequency) . On the other hand, 
when discrete-time systems are analyzed, it is often more convenient 
to express variance in terms of the total mean-square value. The two 
points of view are equivalent provided that the sampling frequency f^ = 

1/T is used for dividing up the power on a two-sided basis. Alternately, 
the folding frequency = f^/2 could be used on a one-sided basis. 

Consider now the process of representing a finite set of samples 
taken with an A/D converter by an ideal train of inpulses. In the 
frequency domain, the spectrum of an impulse train 6^(t) is expressed 
as 


F[«T<:t)] - f E Tifj) C36) 


However, an actual train of nonzero width pulses would have a spectrum 
weighted by d = tj/T, which would approach zero as tj 0. Consequently, 
the impulse spectrum has been effectively weighted by 1/tj by using 
the ideal impulse approximation. On the other hand, when the sampled 
signal is converted back to a continuous- time signal, it can be thought 
of as being held by a zero-order holding circuit whose magnitude is 
unity. The transfer function for this device can be expressed as 
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tj sine fti so that the spectrum is weighted by tj, the net effect 
is tix(l/ti) = 1 , and the level of the signal is preserved. 
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SAMPLING EFFECTS OF NOISE SPECTRUM 


A major objective of this study was to identify possible schemes 
by which the processing and reduction of radiometer data could be 
achieved by digital signal— processing techniques. While the general 
field of digital signal processing is quite mature, there has been 
relatively little work reported in the literature on its application 
to radiometer type signals and systems. There are some peculiar 
aspects of the radiometer process that result in difficulty in 
applying standard signal-processing techniques directly to the signal. 

A significant portion of this study was devoted to investigating the 
nature of the signal and possible ways of circumventing these 
difficulties. 

The "signal" at the front end of a microwave radiometer is 
actually wideband noise having a Gaussian amplitude distribution with 
zero mean. The wideband input filter establishes the bandwidth over 
which the measurement is performed, which typically ranges from 
20 MHz to 2GHz. The square- law detector changes the form of the signal 
to one with a chi-square distribution with one degree of freedom. The 
mean of this signal is proportional to the total power in the input 
signal. Along with the mean value, which represents the quantity to 
be estimated, is a large background noise signal which must be smoothed 
to obtain the value of the true estimate. This undesired background 
signal has a triangular power spectrum. Subsequent processing of 
the detector output reduces significantly the level of the background 
while providing an estimate of the true mean. 

In a Dicke radiometer, correlation of the signal in order to cancel 
out the effective receiver temperature noise is also achieved in the 
postdetector processing. In the case of a noise-injection feedback 
radiometer, an additional objective that must be performed is the 
development of a low-frequency control signal to provide a proper feedback 
mechanism for the loop. 
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One of the earliest questions of interest in investigating 
possible digital signal-processing schemes is to determine an 
appropriate point for sampling the analog signal. At first glance, 
one might conclude that there would be no hope of sampling the analog 
signal prior to the square- law detector since the bandwidths there 
could be as high as 2 GHz and hence result in a required sampling 
rate of 4 GHz using the conventional sampling rate interpretation. On 
the other hand, one might argue that since the mean value, not signal 
reproduction, is the only quantity of real concern, perhaps an under- 
sampled signal could be usable. After all, the Dicke process itself 
as utilized in the front end of the radiometer is a form of sampling 
and its rate is far below the required Nyquist rate. 

As a result of these questions, an analytical investigation was 
made to determine the possible implications of undersampling a wide- 
band noiselike signal. Certain simplifying assumptions were made in 
order to keep the analysis to a manageable level that would provide 
some useful practical interpretations. 

Assume a very wideband noise source n^(t) whose power spectrum 
is G^^(f) = n/2 on a two-sided basis as shown in figure lS(a) . Assume 
first that the signal is filtered directly by a low-pass filter with an 
equivalent one-sided noise bandwidth as shown in (b) . The instan- 

taneous output noise of this filter is noi(t), and the total variance 
associated with this noise is Nqx- From the basic definition of noise 
bandwidth, this variance is simply 

Noi = 7 = 

Now consider the system shown in (c) of figure 18. The signal is first 
sampled (or modulated) by a pulse signal p(t) whose form is shown in 
(d) . Consequently, can be expressed as 

T^sCt) = n^(t) p(t) (38) 


27 


Since p(t) is periodic, it can be expressed as a Fourier series of 
the form 


P(t) 



j 

e 


no) t 
s 


(39) 


where f^ = sampling rate and T is the time between 

successive samples. For the pulse train p(t) with unity pulse ampli- 
tude as shown, the coefficient c^ is given by 

c = d sine nd (40) 

n 

where d = tj/T is the duty cycle of the pulse train. Substitution of 
equation (39) in equation (38) yields 


n (t) = n. (t) 
s 1 



j 

c e 
n 


noj tl 



(41a) 


00 
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n 


n^ (t)e 




(41b) 


The autocorrelation function corresponding to ^g(t) is 


R (t) = E fn (t + y) n*(t + y - t)1 (42) 

ns Ls s -* 

where y is uniformly distributed over the range 0 _< y £ T and is a 
variable introduced as a more rigorous means for establishing stationarity 
of the process. expanded as follows: 
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ns 


(t^) = + w) K 


jnw (t + p) 


c e 
n 


] 


[ ~ -jnw (t + p - t)"I\ 
(t + y - t) 2-< c* e J | 


( 43 ) 


Interchanging the order of the averages results in 


ns 


(t) = [n^(t + y) n^(t + y - t)] 


jnw (t + y) 

Vc e ® 
n 



(44) 


The power spectrum 
by computing the Fourier 


G g(f) of the sampled signal is determined 
transform of equation (44) . 


G (f) 

ns 





G . (f - nf ) 
ni s 


E 



J1 

2 


oo 



(45a) 

(45b) 


(45c) 


where C is the summation. In comparison with the unsampled signal, the 
power spectrum of the sampled signal is multiplied by the infinite 
summation shown. 
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The nature of this summation will now be investigated. It can 
be manipulated as follows: 

C = sinc^nd 

— oo ..oo 





( 46 ) 


The summation in equation (46) converges to a value of tt. Using this 
result, equation (46) may be expressed as 

C = - X 7t = d (47) 

n 

Finally, the output noise power Nq 2 under this condition becomes 

No 2 = d y (2B|^) = dnBj^ (48) 

Comparing equation (48) with equation (37) , it is seen that the 
output noise fluctuation power is d times the value that would be 
obtained without sampling. At first glance, this sounds appealing, 
since the total fluctuation power is apparently reduced by the sampling 
process. However, to truly evaluate this approach, it is necessary to 
consider the existence of a signal component as well. If such a signal 
component is present, it appears as a frequency domain impulse at dc. 
Let Vi represent the value of the signal component (mean value to be 
estimated). The signal power Sqi of the mean value for the continuous 
nonsampled case is simply 

Soi = (Vi)2 (49) 
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With sampling, it is easy to show that the dc component V 2 is simply 

V 2 = dVi (50) 

The signal power Sq 2 with sampling is then 

S 02 = d2(Vi)2 = d^Soi (51) 

The ratio of the signal-to-noise ratio after sampling to the signal-to- 
noise ratio before sampling is determined by 


So2/Nn2 ^ Nqi ^ j2 . i h 
Soi/Noi Sqi Nq 2 d 


(52) 


Thus, the signal-to-noise ratio is degraded by a factor of d when the 
wideband noise signal is undersampled. 


A further explanation of this process will be made. The noise 
components are uncorrelated; therefore, the fluctuation power reduces 
in direction proportion to the time of observation. On the other hand, 
the signal (mean value) is a coherent component and its voltage level 
reduces in direction proportion to the time of observation, meaning that 
its power reduces in proportion to the square of the time of observation. 
Thus, the signal power decreases at a much faster rate than the noise 
power . 


The preceding development has been idealized somewhat in order to 
illuminate the concept fully. The noise was considered separately from 
the signal to simplify the derivation. When both the mean value (signal) 
and fluctuation noise appear together, an additional term appears in 
equation (44) , which is the square of the mean value multiplied times the 
quantity in the brackets. The corresponding spectrum of equation (45c) 
then contains an infinite number of line components at integer multiples 
of the sampling frequency. Depending on the manner in which the result 
is filtered in order to enhance the signal with respect to the noise, these 
line components may or may not affect the net overall fluctuation level. 

The exact ratio given by equation (52) applies only when these line 
components are eliminated. However, the signal-to-noise ratio is degraded 
by the process of viewing the signal for a portion of the time when the 
signal is undersampled. 
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When the sampling rate is sufficiently high that no aliasing 
occurs, only the c^ term in equation (45c) contributes to the noise 
at low frequencies. In this case, the noise power, like the signal 
power, is proportional to and both quantities decrease at the 

same rate. With proper filtering in this case, it is possible to avoid 
any degradation of the signal-to-noise ratio. 

An interesting inference from this development is that the basic 
Dicke radiometer fits the aliasing result as a special case. With the 
Dicke radiometer, the signal is effectively sampled for a half-cycle, 
and thus d = 0.5. Under optimum conditions of estimation, the noise 
fluctuation power-to-signal power has an additional factor of two, 
which agrees with the preceding results. 

From the preceding development, it can be deduced that undersampling 
a noise spectrum is not, in general, an advisable procedure in digital 
processing except when special objectives are sought, such as in the 
Dicke radiometer. 
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OPTIMUM SAMPLING RATE 


Analysis of samp led- data signal fluctuations can be approached 
from. two separate points of view; (1) a frequency domain approach based 
on filtering of the signal and (2) a purely statistical point of view 
based on a reduction of the variance. Both points of view are correct if 
properly interpreted. Furthermore, it appears that the frequency 
domain point of view is most amenable to interpretation when the signal 
is oversampled, while the statistical point of view is more convenient 
when the signal is undersampled. The approaches are about equal in 
complexity when the signal is sampled exactly at the minimum Nyquist 
rate. 

The development that follows should lead to some interesting and 
useful interpretations for these concepts. Consider a signal whose 
two-sided power spectral density G^(f) is of a rectangular form 
as shown in (a) of figure 19. This function is then described by 

Gj^(f) = ^ < f < B. (S3) 

= 0 otherwise 

The total variance of the signal is 

= 2B^ _n = B^n (54) 

Now assume that the signal is sampled by an ideal impulse sampler 
with sampling rate f^. For convenience the impulse sampler will be 
assumed to have a relative weight T = so that the impulse train 

6^(t) can be expressed as 


00 

6„(t) = T Z 6(t - nT) 

1 — oo 


(55) 
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The reason for this assumption is to cancel out the 1/T factor that 
appears as a multiplier for the spectral components of the sampled 
spectrum of an ideal impulse train, which would create an additional 
"confusion factor" in interpreting some of the results that follow. 
This factor actually does not affect the final results anyway, since 
the holding/reconstruction circuit effectively cancels the 1/T factor, 
but by putting in the T factor as assumed, the levels of the 
spectral components can be better interpreted. 

After the signals are sampled, a sum-and-dump algorithm will be 
applied to the set of discrete numbers involved. This algorithm reads 


Y = 


1 N 


n=l 


(56) 


where represents the set of N discrete output values of the sampler, 
and y represents the estimate obtained from application of the algorithm. 
The average is performed over an interval of t seconds, where 


T = NT 


(57) 


Three possibilities will now be considered: (1) Nyquist rate sampling, 

(2) under samp ling, and (3) oversampling. 

Nyquist Rate Sampling 

For this first case, the sampling rate will be assiomed as 

f = 2B. (58) 

SI 

The form of the sampled spectrum Go^^^(f) is shown in (b) of figure 19 
Observe that there is neither aliasing nor spectral gaps. The spectrum 
has exactly filled in the gaps without any "overcrowding." 

From a frequency domain point of view, the effect of the sum-and- 
dump algorithm can be represented by an equivalent noise bandwidth 
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which on a two-sided basis can be shown to be 


B 


N 


1 = J_ 

T NT 


( 59 ) 


The variance 


of the estimate 

y 


y 


can then be expressed as 


nf 

.2 _ Hr = JL = _J]_ = i 

y " 2 N 2t 2NT 2N 


(60) 


The several forms listed provide different interpretations, but they are 
essentially frequency domain in form. These results relate directly to 
the traditional continuous integrate- and- dump circuit in which the output 
variance is determined by a similar form. 

A different interpretation of equation (60) is obtained by noting 

that 



s 



(61) 


Substituting of equation (61) in equation (60) results in 



This form is the familiar statistical result concerned with the sample 
mean concept. The variance of the sample mean is the variance of the 
process sampled divided by the number of samples, provided that the 
samples are all statistically independent. It is seen that this concept 
agrees exactly with the frequency domain point of view when the signal 
is sampled at the Nyquist rate. However, results when the signal is 
either undersampled or oversampled will be considered next. 


Unders amp ling (f^ < 2B^) 

Although the results that will be developed generally apply to any 
degree of undersampling, it is convenient to assume that the bandwidth 
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is an integer multiple of the folding frequency. Thus, let 



( 63 ) 


where M is an ’’undersampling" integer factor. For = 1, this case 
reduces to the Nyquist sampling case. However, the case of interest 
involves > 1. 

Because the signal is undersampled, aliasing of the spectrum will 
occur. The resulting sampled spectrum Go^^^(f) Will be of the form 
shown in figure 20. (This case corresponds to = 4.) The aliasing 
results in exactly components adding together in the low frequency 

range. Since these components add incoherently, the fluctuation power 
in the low-frequency range is exactly times the value of the unsampled 

signal . 

From a frequency domain point of view, the variance is now 


a 


2 


y 


nMy 

2 


nMu _ nMy _ nM^f^ 
2 T 2NT 2N 


(64) 


Comparing equation (64) with equation (60), it initially appears that 
the variance has been increased by a factor over the case of Nyquist 

sampling, and this result is correct if all the other factors in the 
expression are the same as before. However, the input spectral density 
n is in a sense now divided over a greater bandwidth for a given input 
variance. The input variance is 



Substitution of equation (65) in equation (64) results in 


a 


2 


y 


a?M f a? 
L u s _ 1 

2B.N IT 
1 


where the assumption of equation (63) is utilized. 


( 66 ) 
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Hence this result is exactly the same as equation (62) for the 
case of Nyquist sampling because of the fact that all samples are 
statistically independent. Thus, the question of whether the output 
variance is greater or not with undersampling depends on the point of 
view. From a purely statistical point of view, the output variance is 
1/N times the input variance and the reduction effect is the same 
whether Nyquist sampling or undersampling is used for the given number of 
points. However, from a frequency domain point of view, the fluctuations 
appear to be greater; but, for a given input variance, the power density 
is smaller, representing the fact that noise power is spread over a 
wider bandwidth than before. The fact is, however, that undersampling 
is not optimum in the sense that a greater reduction of variance could 
have been achieved in the same time by employing more samples. In other 
words, undersampling results in the variance being reduced as much as can 
be expected for the number of points involved, but it could be reduced 
more in the same amount of time by using more samples . 

Oversampling (f^ > 2B^) 

As in the case of undersampling, an integer relationship between the 
bandwidth and the sampling rate will be conveniently chosen. In this 
case, the following form will be assumed; 


f 

s 


= M (2B. ) 
o^ 1 


(67) 


where is an "oversampling" factor and is the ratio of the actual 

sampling rate to the Nyquist rate. The resulting sampled spectrum will 
be of the form shown in figure 21. (This case corresponds to = 2). 
Note in this case that there are "holes" in the resulting spectrum. 

The output variance in this case is 


nf 

.2 _ Hr = n _ n _ s 

y 2 N 2 t 2NT 2N 


( 68 ) 
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A statistical intepretation of this result is obtained by substi- 
tuting f^ from equation (67) in equation (68) and utilizing the fact 
that 



The result is 


a 


2 

y 



(70) 


In this case, it appears that the output variance is increased by a 
factor of over that of the Nyquist sampled case, and this point 

of view would be correct if a? were the same as for the Nyquist 
sampled case. However, the input fluctuations are nov\r confined to a 
narrower bandwidth for a given input power spectrum. From the fre- 
quency domain point of view, the variance has been reduced as much as 
can be expected for the given integration time. In this case, the 
frequency domain point of view seems to be easier to interpret due to 
the lack of aliasing in the spectrum. 

To summarize this case, oversampling results in the variance being 
reduced as much as can be expected for the particular time involved, but 
it could be reduced more by letting the same number of points extend 
over a longer period of time . In other words, the samples are too close 
together to be completely statistically independent, and increasing the 
time between samples will result in a greater reduction of the variance 
for the same number of points. 

From the preceding developments, it can be concluded that the 
Nyquist rate is an optimum rate for the reduction of variance in a sampled- 
data signal with rectangular power spectrum both from the frequency domain 
point of view and from the statistical point of view. At that rate, 
both the filtering and the statistical independency properties are "doing 
their best." 
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The optimum rate for this purpose should not be confused with 
the signal reconstruction objective in other applications where a rate 
exceeding the Nyquist rate is always used. As a matter of fact, the 
assumption of a flat spectrum has been made thus far, so it is logical 
to ask what modification is necessary when the power spectrum rolls 
off more gradually, as it typically does. 

Consider, for example, the gradual spectral rolloff shown in (aj of 

figure 22. In order to make the aliasing error vanishingly small, the 

sampling rate f is selected to be as shown in (bj . While the 
Si 

resulting increase in the aliasing power contained up to the folding 
frequency is negligible, the variance reduction is suboptimal in that 
it is less than the 1/N factor. Specifically, the output variance 
resulting from a sample mean definition can be most directly calculated 
from a frequency domain point of view and is 



(71) 


An alternate strategy with such a rolloff as this is to undersample 
somewhat so as to ’’conserve the sampling rate" with the philosophy that 
the variance won't increase much in the process. The results of this 
assumption are shown in (c) of figure 22. The output variance in this 
case is 


r2 _ 


0.5£ 

/ 

*'-0.5f 



(72) 


where G 2 (f) is the modified spectrum with aliasing. 
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The case where aliasing is present (fig. 22c) results in an 
increase in the output fluctuations as compared to the case of figure 
22(b), but the decrease is more "efficient" for the same number of 
points. The tradeoff of interest is whether or not the reduction in 
the required sampling rate is worth the cost of the increased variance. 
Each case would have to be considered on its own merits. If the 
data rate is rather slow and the digital circuitry is capable of 
operating at a high rate anyway, the optimum choice will likely be to 
sample at a high enough rate to eliminate all aliasing error. 
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SECOND-ORDER DIGITAL LOOP DESIGNS 


Early in the course of this study, several second-order control 
loop designs were investigated and successfully simulated. The choice 
of a second-order loop was made at that time because the current analog 
system utilized a second-order loop. As will be shown later, a first- 
order loop followed by post- loop sum-and-dump filtering was finally 
selected for the prototype implementation. Nevertheless, there are 
still some worthwhile features of the second-order loop that could make 
it a viable candidate for some applications. Consequently, some of 
the intrinsic features of a second-order digital loop will be 
discussed in this section, and some representative designs will be 
shown. 

The concept of utilizing the second-order loop was based on 
performing virtually all the data smoothing in the loop itself and con- 
sidering little or no post-loop filtering. It was also decided to employ 
coefficients in this design that could be realized exactly by a rather 
limited number of bits. In this case exact coefficients were known, 
making the equivalent smoothing bandwidths more readily predicted. 

In establishing equivalence, or at least correspondence, between 
a digital loop and an analog loop, several approaches are possible. 

None of these are "exact” in one sense since the discrete-time system 
is different from the continuous -time system, and the best that can be 
done is to establish a correspondence with respect to some particular 
criterion. For that reason, a comparison with respect to more than one 
criterion will be made in some of the developments that follow. 

One criterion is the direct z- transform correspondence z = or, 
equivalently, s = (1/T) £n z. This relationship computes the values of 
s and z that correspond to each other through the z-transform 
definition, and in this sense is exact, 

A second criterion is a correspondence through the bilinear 
transformation. A popular and convenient means for designing digital 
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filters is to employ a prototype continuous system variable p, obtain 
a continuous transfer function approximating the desired behavior, and 
then to replace p by 


P = 


K(1 - z-^) 
1 + z“^ 


(73) 


where K is a mapping constant. For very low frequency correspondence 
between the prototype analog frequency and the final digital frequency, 
it can be shown that an optimum choice of K is 

K = I = 2f^ (74) 


The actual second-order functions chosen for the loop were based on 
a modified bilinear transformation criterion. Consider momentarily a 
first-order analog function of the form 


G(p) 


a 

p + a 


(75) 


When the bilinear transformation is applied to equation (75), the 
corresponding H(z) is determined as 


H^(z) 


a 

K + a 


(1 + z~\) 
1 - Bz"^ 


(76) 


where 


K - a 

K + a (77) 


Consequently, if B is selected as a convenient value in the digital 
filter, the corresponding a in the analog sense can be determined as 


K(1 - B) 
1 + 


( 78 ) 
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The first design of a digital estimating circuit considered is 
shovfn in figure 23. The particular structure shown in the figure 
represents only the digital processor portion starting from the A/D 
converter and continuing to the point at which the output estimate is 
determined. The feedback noise-injection scheme is not shown in the 
figure. The difference equations describing this system are as follows: 

u(n) = x(n) + x(n - 1) + u(n - 1) (79) 

v(n) = 2"^ u(n) (80) 

w(n) = v(n) + v(n - 1) + 0.875 w(n - 1) (81) 


y(n) = w(n) 


(82) 


The corresponding transfer functions are 


U(z) _ 1 + z-1 
X(z) = 1 _ ,-l 


(83) 


V(0 _ 2-3 
U(z) ■ 


(84) 


W(z) ^ 1 + z 

1 - 0.875z~l 


(85) 


y(z) _ 2-3 

W(z) 


( 86 ) 


The composite transfer function is 


„f_l Y(z) _ 2-6(1 . z-l)2 

^ (1 - z-l)(l - 0.875Z-1) 


(87) 
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The low-frequency behavior of HCz) may be determined by substi- 
sT 

tuting z = e and letting s become very small. The function H(e ) 
then asymptotically approaches 


H(e 



2 ^(2)2 _ 0.5 
sT(0.125) sT 


(for small s) 


( 88 ) 


For this particular design, with the sampling rate selected as f^ = 750 Hz, 
the function approximates 375/s. 

Based on the constant B = 0.875 in equation (76), the corresponding 

CtT 

analog a may be approximated by either setting e~ = 0.875 or by 
using the bilinear transformation. In the first approach, a = 100.15 
rad/s. In the second approach, a = 100 rad/s based on equation (78). 

The two values are in close agreement, and correspond very closely to 
the optimum analog value selected in an earlier redesign effort. 

The actual analog transfer function G(s) of the estimation circuit 
was earlier selected as 


rr = 17.857 X 8.82 
“ s(l + 0.01s) 


For small values of s, this function approximates 


(89) 


G(s) « iilil (90) 

s 

The corresponding digital function has more than twice the gain required, 
so the signal should be attenuated by a factor of 157.5/375 = 0.42 before 
processing. 

The second approach considered in second-order loop designs was one 
in which two separate modes could be employed. The first mode, designated 
mode A, is the lower resolution, faster response mode. On the other 


44 



hand, mode B provides a greater resolution at the expense of a 
longer integration time. For this particular design, an Intel 8086 
microprocessor development unit had been identified for prototype 
development work. Consequently, its clock frequency of 4.9 MHz was 
selected for timing, and the sampling rate was chosen as 


f = 4.9 X 10^/213 = 598.14 Hz 
s 

The discrete transfer function of the lag- lead digital filter was 
selected in the form 


H,Cz) 



(91) 


The value of B for mode A was selected as 

B = O.IIOI 2 = 0.8125io (92) 

where the subscripts 2 and 10 represent the binary and decimal bases, 
respectively. The corresponding value of a in the analog representa- 
tion is determined from the direct z-transform definition as 

a = -(1/T) £n 0.8125 = 124.20 rad/s 

and from the bilinear transformation as 

a = 2 X 598.14 x(l - 0.8125)/(1 + 0.8125) = 123.75 rad/s 

To simplify further references, the value a = 124 rad/s for mode A will 
be assumed in subsequent discussions. 

The layout of mode A is shown in figure 24. The various difference 
equations are tabulated as follows: 
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u(n) = x(n) + x(n - 1 ) + u^n - 1 ) 


(93) 


v(n) = 2 "‘^u(n) 


w(ti) = v(n) + v(n - 1) + 0.8125 w(n - 1 ) 


y(n) = 2 "^ w(n) 


The corresponding transfer functions are 


U(z) _ 1 + z"l 
X(z) = 1 _ ,-l 


V(Z ) _ .-4 
U(z) 


W(z) ^ 1 + z-^ 

1 - 0.8125z"l 


IM. - 2-2 

W(z) 


The composite transfer function is 


(94) 

(95) 

(96) 

(97) 

(98) 

(99) 

( 100 ) 


Y(z) ^ 2 ~ 6 (i + z-l )2 

(1 - z“l)(l - 0.8125Z-1) 


( 101 ) 


As observed in equation (90) , the analog estimating circuit approxi- 
mates 157.5/s for small s. The function of equation (101) approximates 
2" ^(2) 2 /sT( 0. 1875) = 199.38/s. Thus, a small attenuation of the analog 
signal before sampling should result in nearly equal behavior of the 
digital system as compared with the analog system. 
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The value of B for mode B was selected as 


B = O.IIIIOI 2 = 0.953125io (102) 

The corresponding value of a from the direct z-transform was computed 
as a = -(1/T) £n 0.953125 = -28.716 rad/s and from the bilinear transform- 
ation as a = 2 X 598.14(1 - 0.953125)/(1 + 0.953125) = 28.71. The 
latter value will be selected for subsequent references. 

The layout for mode B is shown in figure 25. The various difference 
equations are tabulated as follows: 


u(n) = x(n) + x(n - 1 ) + u(n - 1 ) 

(103) 

v(n) = 2 ~^x(n) 

(104) 

w(n) = v(n) + v^n - 1) + 0.953125 w(n - 1 ) 

(105) 

y(n) = 2 “^ w(n) 

(106) 

corresponding transfer functions are 


U(z) 1 + z-l 

(107) 

V(z) _ 2-5 
X(z) 

(108) 

W(z) _ 1 + 

1 - 0.953125Z"! 

(109) 

Y(z) _ 5 

W(z) - 

( 110 ) 
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The composite transfer function is 


^ (1 - z"^)(l - 0.953125z“l) 

The function of equation (111) approximates 2”^^ (2)^/sT(0. 046875) = 
49.845/s for small s. The corresponding analog filter transfer function 
was quite close when redesigned for this higher accuracy mode. Note 
that the present analog system does not operate in a higher resolution 
mode such as this. Instead, higher resolution is achieved by post- 
loop processing of the data obtained from the lower resolution operating 
mode. 
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CANCELLATION OF DICKE RIPPLE 


A troublesome problem in an analog Dicke radiometer system is 
the presence of the Dicke ripple component in the data output. 

Because of the periodic switching in such a radiometer, an undesired 
disturbance at the switching frequency appears in the signal output. 
This disturbance consists of a fundamental component at the switching 
frequency plus components at odd integer multiples of the switching 
frequency if symmetrical switching is employed. Due to the normal 
low-pass nature of the forward transfer function, tlie most troublesome 
component is the fundamental, and its magnitude is a close approxima- 
tion to the total ripple level in most cases. 

While it is possible to set the loop parameters to adjust this 
ripple to a tolerable level by careful design, an interesting concept 
came to light in investigating possible digital-processing schemes. 
Consider the discrete transfer function H(z) of the trapezoidal 
integration approximation as given by 


H(z) 


K (1 + z ~^3 
1 - z~ ^ 


C112) 


The steady-state transfer function corresponding to this function 

is readily shown to be 




ir (J^T 

K cos 
sin -x- 


= K cot nfT 


(113) 


The form of the magnitude response [H(e^^^)| corresponding to equation 
(113) is shown in figure 26. The frequency response approximates that 
of an ideal integrator in the very low frequency range, i.e., for f << 1/2T. 
However, for f = 1/2T, the steady-state transmission is identically 
zero. 
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Assume now that the sampling frequency - 1/T is chosen to be 
twice the Dicke frequency f^, i.e., f^ = 2f^. This results in a 
zero or null of transmission at the Dicke frequency and at odd 
integer multiples of the Dicke frequency. Thus, as long as the 
switched signal satisfies the half-wave symmetry conditions of Fourier 
theory, all undesired components in the ripple disturbance will be 
completely eliminated. 

As it turns out, the proposed feedback schemes utilize a pulse- 
injection process which disturbs the half-wave symmetry of the switched 
signal. The resulting signal in this case does contain even harmonics 
which are not eliminated by this process. However, such components are 
easier to handle since they are higher in frequency and smaller in 
magnitude. The major point of interest is that a significant reduction, 
or possibly an elimination, of the switching ripple can be achieved 
directly through a digital algorithm. 
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PROPOSED CONTROL LOOP 


After investigating in some detail the several second-order 
loop designs previously discussed, attention was directed toward 
the concept of a first-order loop followed by post- loop sum-and-dump 
filtering. Because of the presence of the post-loop filtering, the 
system was actually of higher order, but for subsequent references the 
term "first-order processor" or simply "first-order loop" will be 
used. 

Some of the advantages of the first-order loop determined from 
this investigation are (1) the form of the response is less sensitive 
to parameter deviations than the second-order loop. Small variations 
in gain parameters can sometimes adversely change the damping ratio of 
a second-order loop. C2) The loop can be made to respond faster, 
thereby allowing a greater reduction of fluctuations by the sum-and-dump 
post-loop filter. (3) It is easier to take advantage of the Dicke 
cancellation scheme previously discussed when the loop is a first-order 
design. 

The basic design of the first-order estimation loop filter is 
shown in figure 27. The difference equations for this system are 


w(n) = x(n) + x(n - 1) + w(n - 1) (114) 

y(n) = 0.25 w(n) (115) 

The corresponding transfer function H(z) for the filter is 


_ Y(z) _ 0.25(1 + z~l) 
- X(z) - j ^ 


C116) 


The design of the loop estimation and feedback system is shown in 
figure 28. The timing references for this system are derived from the 
4.9 MHz and the 2.45 MHz signals available in the microprocessor unit. 
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The sampling rate for the analog signal is = 2.45 MHz/2^^ = 598.14 Hz. 
The A/D converter is an Analog Devices type AD-572 12-bit successive 
approximation unit. The drive for the Dicke frequency f^ is 
obtained by dividing f^ by 2, thus giving f^ = 299.07 Hz. This 
results in zeros of transmission at odd integer multiples of the Dicke 
frequency as previously discussed. 

The analog aliasing filter is a cascade of 3 simple one-pole 
filters, each having a 3-dB break frequency at 100 Hz. An investigation 
was made to determine the need for more sophisticated aliasing filter 
designs having complex poles, inasmuch as this particular filter has a 
rather pronounced amplitude ’’droop" in its passband. However, since the 
only component of primary interest is the dc component, and since the 
actual passband shape is relatively unimportant, this filter was found 
to be perfectly adequate, and it simplifies the overall design. There is 
a small aliasing error, but its level is insignificant in view of 
other error contributions. 

The noise-injection scheme proposed for this design differs 
considerably from the original design in that it will be injected 
continuously during a portion of the Dicke cycle. In this manner, the 
bandwidth requirements of the broadband amplifier preceding the Dicke 
correlator are eased considerably. The noise is injected once per 
cycle as illustrated by figure 29. (Actually, it is injected on both 
half-cycles, but since the Dicke switch only sees it once per cycle, 
the extra injection does not appear on the figure). The length of time 
that the noise-injection pulse is on is determined by the value of the 
digital word at the output of the loop digital filter. A word is 
loaded in the down-counter once per sample time T. The noise diode is 
gated on through the PIN switch at the beginning of this interval. When 
the down-counter reaches zero, the "status" signal turns off the switch, 
thereby reducing the noise injection to zero. 

The maximum duty cycle corresponds to minimum input temperature, 
which results in maximum feedback temperature injection. Based on a 
reference temperature of 308 K, a directional coupler coefficient of 
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-20 dB, an excess noise temperature of 365,000 K, an additional attenua- 
tion of 10 dB in the feedback loop, and an input temperature of 0 K, 
the maximum duty cycle was found to be 0.835. Conversely, the minimum 
duty cycle corresponding to 300 K input was found to be 0.0217. The 
range of temperatures and duty cycle are illustrated on figure 30. 
Operation near minimum duty cycle poses the same sort of problems as 
the current analog system since the resolution is degraded seriously 
in that range. 

It was decided not to operate at a duty cycle of unity due to 
"spillover" resulting from finite bandwidth limitations. Consequently, 
the duty cycle of 0.835 was judged to be a good choice for the 
maximum value of the duty cycle. 

While odd harmonics of the Dicke frequency are eliminated by the 
digital filter, even harmonics will appear due to the asymmetrical noise 
injection. Referring to figure 29, an analysis of the levels of even 
harmonics present may be made from a combination of the energy balance 
and Fourier theory. At a condition of zero steady-state error, the 
area of the positive half-cycle must equal in magnitude the area of the 
negative half-cycle. This means that 

C Tjti + C(T^ + T^)T = C(Tg + Tj^jT (117) 

which results in 

Tjti = (Tg - T^)T ( 118 ) 


or 


Tjd = Tg - T^ (119) 

where d = ti/T is a duty cycle defined as the ratio of the "on" interval 
to half of the Dicke period. (Equivalently, it is the total on time. 
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including the superfluous extra injection, per Dicke period to the 
total period.)* Full injection corresponds to d = 1. 

The even harmonics that arise from the asymmetry may be considered 
to originate from a pulse train having a period 2T, a duty cycle d/2, 
and an amplitude CTj. The magnitudes of these components on a 

conventional sine-cosine basis can be expressed as 


A 

n 



sin mrd/2 
n7rd/2 


= C(Tg 




sin 


nTTd/2 

md/2 


( 120 ) 


Since the fundamental and odd harmonics are eliminated by the loop 
filter, the most significant component of concern is the second harmonic 
term. This value is 


^ I IT 

As one would expect, A 2 = 0 for d = 1. Conversely, the maximum value 
of A 2 occurs for d = 0.5, meaning that noise is injected half of the 
time between successive analog samples. This maximum value A 2 = CT^/tt 
occurs for an input temperature ^ 124 K. This could be serious 
except for the fact that the 3-pole aliasing filter has an attenu- 
ation of about 47 dB at the sampling frequency, plus the rolloff 
of the loop itself provides an attenuation of about 37.5 dB. The overall 
level of the even harmonics was found to be negligible in all the 
simulations . 

A simplified block diagram of the control model for the discrete- 
time system is shown in (a) of figure 31. The gain of 4.61 x 10 ^ 
represents the net effect of the original gain of 2.88 x 10 ^ times an 
additional gain of 1.6 that was later added to establish an optimum 
transient response. The second block is a "redefinition" gain 
established by the following procedure: The input range of the A/D 

converter is set at ±2.5 V, which was determined by simulation to be 
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the proper level. This means that the least significant bit at the 
input corresponds to 5V/2^^ = 1.22 mV. In passing through the 
processor, a 4 to 1 expansion of the signal level occurs, meaning that 
two additional bits are required. However, it is desired to "redefine” 
the signal on a normalized basis so that each digital word has a 
maximum value of unity. Considering the 2.5 V maximum at the input 
coupled with the 4 to 1 expansion, the output would really have a 
maximum of 10 V on an analog basis, so a "gain" factor of 0.1 is used 
to establish it at unity maximum level. 

The forward gain function Hj (z) is then 

Hi(z3 = 115.25 X 10"^ ^ (122) 

1 - z“ ^ 

For very small s, this function behaves approximately like an analog 
gain function 


Gi(s) « 


115.25 X 10"®(2) 
sT 


0.158 

s 


(123) 


Tiiis equivalence is shown in (b) of figure 31, which clearly 
illustrates the nature of the corresponding first-order analog system. 

The appropriate loop gain content is 0.138 x 365 = 50.32, which will 
be rounded off to 50 for subsequent calculations. The closed loop transfer 
function G(s) will then be expressed as 


G(s) 


0. 157 ^ 2.74 X IQ-3 
s + 50 1 + 

50 


The overall proposed system is shown in figure 32. 
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POST- LOOP PROCESSOR 


The digital filter contained in the loop was designed primarily to 
provide the necessary dynamic loop behavior for optimum control of 
the noise-injection process. In contrast, the primary objective of 
the post- loop filter was to provide a smoothed estimate of the bright- 
ness temperature in as short a time as possible commensurate with the 
sensitivity specifications. 

The concept used to reduce the variance of the temperature estimate 
in the post-loop processor was the well-known sample mean algorithm. 

This process will hereafter be denoted as a "sum-and-dump" algorithm 
due to its close similarity to the integrate-and-dijmp circuit used in 
analog matched filter and estimation systems. Indeed, the mathematical 
behavior of the sum-and-dump algorithm on a discrete-time basis is 
virtually identical to the behavior of the integrate- and- dump filter 
on a continuous-t ime basis. 

The integrate-and-dump analog filter is optimum in that the 
product of the equivalent noise bandwidth times the settling time is 
minimum for all analog lumped filters. The corresponding sum-and-dump 
algorithm possesses the same optimum property among the class of discrete- 
time filters. 

The post- loop processing is actually a set of several similar sum- 
and-dump algorithms of different lengths as shown in figure 33. Each 
algorithm involves a scanning window that picks up all applicable loop 
samples over the duration of the loop. The number of points chosen in 
a window is always selected as an integer power of two. This simplifies 
the multiplication and division required by the microprocessor, and it 
simplifies the structure of the dependency between shorter and longer 
versions of the algorithm. 

Let y(n) represent the discrete values appearing at the output 
of the loop, and let represent the smoothed discrete values 

appearing at the output of different length post-loop algorithms. The 
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integer N represents the number of points in the particular algorithm. 
Note that n is an integer changing at the loop sampling rate while 
i changes only once for every 16 values of n. Thus, 

n = 16i (125) 

or 

i = integer part of (126) 


The first value is simply yi(i), which is not actually a 
smoothed value but is the output of the loop as held for an interval of 
16 periods at the sampling rate, i.e.,T^ = 26.75 ms. Thus, 


7l (i) 


y(i) for i 



(127) 


The remaining values at the outputs of different stages can all be 
expressed directly in terms of yi(i), or they can be expressed in 
terms of smoothed estimates of lower order- All estimates are updated 
every 16 loop samples even though the statistical dependency of 
subsequent samples increases with the order of the estimate. 

A tabulation of the various estimates follows: 


yi(i) + Xl(i - 1) 

X2(i) = 2 

yi(i) + yi(i - 1) + yi(i - 2) + yi(i - 3) 
y 4 (i) = ^ (129a) 

y2(i) + y2(i - 2) 

= 5 (129b) 
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(130a) 


7 

ysCi) = |- yi(i - I) 
1=0 


y4(i) + y4(i - 4) 

= 2 (130b) 

15 

yi6(i) = yiCi - I) (131a) 

1=0 

ys(i) + ys(i - 8) 

= 2 (131b) 

51 

y32(i) =-^2 (132a) 

1=0 

yi6 Ci) + yi6 (i " 16) 

= 2 (132b) 
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y64(i3 = (133a) 

1=0 

y32(i) + y32(i - 32) 

= 2 (133b) 

Additional representations of the sum-and-dump algorithms from a 
digital filter perspective and a z-transform point of view are shown 
in figures 34 and 35, respectively. 

The variance reduction of the sum-and-dump algorithm will now be 
analyzed. First, the assumption of a uniform power spectrum will be 
made, and this will be modified later to include loop frequency dependency 
Consider then the general form of the algorithm as given by 
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C134) 


N-1 
1=0 

Let represent the z-transform of XjsjCi) Yi(z) represent 
the z-transform of yi(n). Transformation of both sides of equation (134) 
yields 


VnU) 


N-1 

1 V' 

N ^ 

1=0 


--I 


YlCz) 


C135) 


The transfer function of the post-loop sum-and-dump algorithm 

can be expressed as 


N-1 


H (z) - 

N*- ■' ■ YiCz) N ^ 
1 = 0 


-I 


(136) 


The steady-state frequency response will be denoted as 

i coT 

o • 

this quantity is determined by substituting z = e in equation (136) . 
(Note that T is the post- loop sample time, not the loop sample time, 
i.e., T^ = 16T.) The result initially is 

N-1 . „ 

^ 2 ^^ ^137) 

1 = 0 


By means of the summation formula 


N-1 

E 


0 


I 


a 



(138) 
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and the exponental definition of the sine function, equation C137) , can 
be expressed as 


H-Cf) 


sin(NwT^2) 
N sin ((joT^/2) 



(139) 


The equivalent one-sided noise bandwidth can be expressed as 




0.5f 


os 


B., = I |H(f)|^df 

0 


(140a) 


=/ 


^■^^os / sin NttT^T 


N sin ttT f 
o 


df 


(140b) 


where f^^ is the post-loop sampling rate, i.e. f^^ = 1/T^ = f^/16. 

For integration purposes, a change of variables was made by setting 
X = T^f j and the integral then becomes 


®N = 


= f/ 

r\ 


0.5 / . „ '2 
sin Nttx 

N sin TTX 


dx 


(141) 


It can be shown that the value of this integral is given by 


B 


N 



-[log2 N + 1] 


(142a) 


1 _ 0.5 

2NT NT 
o o 


(142b) 


60 



Let T = NT^ = total time interval for averaging. Substituting this 
value of T in equation (142) , the one-sided equivalent noise band- 
width is simply 


B 


N 


0.5 

T 


(143) 


This result is exactly the same as for the continuous -time integrate- 
and-dump filter with t as the integration time. Thus, the sum-and- 
dump algorithm for a discrete time signal functions exactly the same 
as the integrate- and- dump filter for a continuous -time signal provided 
that the summation interval in the discrete-time case is equal to the in- 
tegration interval in the continuous-time case. A further subtle assumption 
is that the power spectra for the two cases would have to be compared 
on the same basis. This would imply optimum sampling at the Nyquist 
rate for the discrete-time system as previously discussed. 


Assume momentarily that the input yi(i) has a uniform power 
spectral density (f) = n on a one-sided basis over the frequency 


range from dc to 0.5 f 


so 


The output variance a ^ is then given by 



0 . 5n _ 0 . 5n 

NT ~ T 
o 


(144) 


The fact is, however, that the input noise power spectrum in the 
case of interest is not flat but has the shape of the power transfer 
function of the loop. The loop relative amplitude-squared response 
function A^(f) expressed as a function of steady-state frequency is 
dominated by a single pole and is of the form 


A^(f) 



(145) 
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where fj is the 3-dB loop frequency. As a close estimate of the 
anticipated design value, the damping constant of the loop is about 
a = 50 corresponding to fi = 50/2ir = 7.96 Hz, which was rounded to 
8 Hz, giving 


A2(f) = 



(146) 


An equivalent noise bandwidth based on the combination of the 
loop response given by equation (146) and the sum-and-dump algorithm can 
now be derived. The form for the equivalent noise bandwidth is 


^NO 



sin NttT f \ 

— 2 _ 

N sinirT f/ 
o 


0 


2 



(147) 


This function was modified for integration by substituting x = T^f as 
before. After some simplification, the result becomes 



(148) 


This integral was evaluated numerically for values of N ranging from 2 
to 64 (in integer powers of 2). It is convenient to tabulate the 


integrals in the form 
plus the case for N = 

N 1 

y = 0.3362 



B^,^NT = y. Thes 
NO 0 

2 

4 

8 

0.3459 

0.4177 

0.4586 


e data are shown below 


16 

32 

64 

0.4793 

0.4896 

0.4948 
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Observe that as the number of points increases, y approaches 0.5 or 
®NO ^PP^oaches 0.5 /t, which is the value without the additional loop 
filtering. In this limiting range, the post- loop algorithm reduces 
the spectrum at such a low frequency range that the effect of the 
loop filter is negligible. 
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SENSITIVITY OF PROPOSED DESIGN 


In this section, an analysis o£ the sensitivity o£ the proposed 
loop design will be performed. Although the loop is a hybrid system in 
that it is partly a continuous -time system and partly a discrete-time 
system, the input-output relationships on a continuous -time basis will 
be maintained since it is from this point of view that the results must 
be interpreted and applied. 

As previously shown, the transfer function of the control loop is 
dominated by a single pole and is approximately equivalent to a 
continuous-time system of the form 


G(s)= ^ 

B A 


0. 157 
s + 50 


(149) 


where Y is the output loop estimate before post-loop filtering. The 
steady-state transfer function G(jw) corresponding to equation (149) 
can be expressed as 


G(jo)) 


0. 157 
50 + j w 


2.74 X 10"3 


1 + 


^50 


2.74 X io"3 


1 + j- 


(150) 


where the 3-dB break frequency is rounded slightly to 8 Hz for convenience. 

The digital words from the output of the loop are then processed 
through the post- loop filter. The equivalent noise bandwidth B^^^ for 
the whole system is determined from the results of equation (148) , which 
were tabulated in the last section. 

The sensitivity AT of the closed- loop noise-injection feedback 
radiometer is determined from the equation 


AT = 2(Tg + Tj^) 



(151) 
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where Tg is the reference temperature (308 K for this system) , 
is the input noise temperature of the receiver, is the equiva- 
lent noise bandwidth, and is the input statistical bandwidth. 


For this system B 


NO’ 


can be expressed as 


Y Y 

®NO ^ T ^ NT* 
o 


(152) 


where y is the constant corresponding to the summation interval as 
given in the previous section. 

Substitution of equation (152) in equation (151) along with present 
radiometer values yields 


y BgjNT " 16.1696 x iq3 


where the value T = 26.7496 ms was used. For different values of N, 
o ’ 

the results of equation (153) are summarized below. 

N = 2 



AT = 

6.725 X 

103 

(154) 


V®SI 


N = 4 

AT = 

5.225 X 

1Q3 

(155) 




N = 8 

AT = 

3.871 X 

1Q3 

(156) 


Jb 

^ SI 
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N = 16 


^ 2.799 X 1q3 

V^i 


N = 52 

AT 


2 X 1q3 



(157) 


(158) 


N = 64 


AT 


1.422 



(159) 


The AT values for different bandwidths were computed and are listed 
in table 1. Along with the values of AT, it is also necessary to 
determine the length of time required for a measurement. The total 
measurement time can be represented as the sum of two components 

of the form; 


T, = T + T (160) 

M s 

The quantity t represents the time required to perform the post- loop 
algorithm and is simply t = NT^ as previously discussed. On the other 
hand, is the settling time of the loop, and its value depends on 

the exact criteria of final value closeness chosen for the loop. 

The following arbitrary but reasonable criterion was employed for 
defining the settling time of the loop. The buildup of the output 
estimate of the loop on a continuous -time basis would be of the form 


y(t) = CiCTj - T^) Cl - (161) 

where is a constant. The error Ay based on a finite settling 

time is 


Ay = Ci(Tg - T^)e-^0^ 


(162) 
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A worst-case value of T,. - T. = 308 was selected. The time x was 

BA s 

then calculated such that the error was no greater than AT/10; i.e., 

the settling time error could be no greater than 10 percent of the 

fluctuation error. This criterion results in the following value for 

T as a function of AT. 
s 


T 

S 



3080 

AT 


(163) 


These data are also tabulated in table 1 along with the calculated 
values of the total measurement time. 
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SIMULATION OF PROPOSED DESIGN 

The proposed system was simulated on the digital computer using 
Advanced Continuous System Language (ACSL) and some special programs 
that were developed earlier under this and a previous contract. The 
actual programs have been described in some earlier publications (refs. 

1 and 2) . 

Because it would consume far too much computer time to simulate 
the radiometer with the actual bandwidths employed in the real-life 
system, a decision was made use a simulated input statistical bandwidth 
of 1/(2)^ = 1/256 of the smallest RF bandwidth. The resulting simula- 
tion bandwidth is then 20 MHz/256 = 78.125 kHz. The reason for this 
seemingly odd choice will now be explained. Since the sensitivity of 
the measurement varies as 1/^ ^SI ’ simulation is 

■^256 = 16 times the corresponding value of the actual design at the 
smallest bandwidth, A factor of 16 corresponds to 4 bits in a binary 
representation. While the actual system will employ 12-bit words, it 
was desired in the simulation to have quantization effects appear in 
the same relative level as they would in the final design. This could 
then be readily achieved by employing an 8-bit A/D converter and a 256- 
step down-counter for the noise-injection counter in the simulation. 
Thus, in the simulation, the relative quantization error as compared 
to the AT level would be in the same proportion as for the final 
design . 

Even with the scaling employed, the computer simulations were quite 
expensive, and only a limited number were used. Nevertherless , the 
results obtained were quite good, and they verify the predicted behavior 
quite well. 

The results of four separate runs with all post-loop outputs are 
tabulated in tables 2 through 7. In interpreting these results, bear 
in mind that the limited number of runs coupled with the dynamic nature 
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o£ the statistical parameters result in some high statistical variations 
from run to run. However, the general trends and averages here are 
most significant in establishing the validity of the concept. 

In each of the tables, the following parameters are tabulated: 


"time" = value of time at which averages are computed, 
temperature estimates 1-4 = the tabulated results obtained at the 

times listed. The input temperature is T = 100 K. The output 
is scaled so as to read 0.99(Tg - T^) = 205.92 K , which would 
be the ideal value. 

= estimate of mean value of process obtained by an ensemble average 
for the four runs 

s = estimate of standard deviation of process 

y N 

a = theoretical predicted value of standard deviation obtained 
^ N 

using results of preceding section 

A 

y^, = mean value of estimated means 
N 

^ = mean value of estimated standard deviations 


Table 2 provides a tabulation of the various data for y 2 (n) , 
i.e., the estimate based on two post-loop samples. Note that the mean 
of the standard deviations is 23.99 K, which compares favorably with 
the theoretical value of 24.06 K. 


As N increases, the statistical dependency between averages per- 
formed at different values of time increases so that estimations of the 
values down the column decrease in effectiveness. In other words, 
the greater the resolution of the estimate itself, the less one is 
able to assess the accuracy of the resolution. 

Table 3 provides a tabulation of the results for y 4 (n). The mean 
of the standard deviations is 17.75 K compared with a predicted value 
of 18.69 K. The values given in tables 4 through 7 provide similar 
results for ya(n), yieCn), y 32 (n), and respectively. 

The results are illustrated graphically in figure 36. 
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APPENDIX 


ANALYSIS OF CLOSED LOOP NETWORK 
AND NOTCH FILTER 


by 

John M. Jeffords 
Old Dominion University 
Electrical Engineering Department 



INTRODUCTION 


Old Dominion University and NASA/ Langley Research Center are 
involved in an effort to develop design techniques for improved micro- 
wave radiometers. This report documents part of the analysis 
performed on circuits used to model a portion of the radiometer system 
in support of the continuing design effort. 

A microwave radiometer is a receiver used to estimate the mean 
value of a noise signal in the presence of background thermal noise. 
Although microwave radiometers exist in a variety of configurations with 
different complexities, sensitivities, and accuracies, the current 
interest of this effort is the noise-injection feedback system. In 
this configuration, the output of a controlled noise source is added 
to the input noise and the sum is compared to a stable reference 
noise. Feedback is used to adjust the controlled noise source so that 
a balance is achieved. The input noise is then the difference between 
the known reference noise and the output of the controlled noise source. 

Although digital signal processing is contemplated for the 
microwave radiometers to be developed, a better understanding of the 
current analog systems was desired. The current system evolved through 
a combination of both analytical and empirical techniques. An analysis 
of this system provides a baseline upon which future design can be 
established. 

Although the actual microwave radiometer involves the interaction 
of a number of random processes, the control mechanism of the feedback 
loop can be modeled by a simpler "deterministic" form which approximates 
the loop operation at very low frequencies. 

STATEMENT OF PROBLEM 
The specific objectives of this study were to: 

1. Determine the optimum value of damping factor for the loop to 

minimize the product of noise bandwidth and settling time. 
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2. Extend the concept of settling time as usually given in books on 
control theory to include a more general criterion for precision- 
bounded measurements. 

3. Determine the relative sensitivity of the loop to changes in the 
open- loop gain factor, particularly in regard to the error. 

4. Investigate the use of a lead time constant as it affects error 
and sensitivity of the loop. 

5. Evaluate the performance of a notch filter. 

GENERAL APPROACH 
Loop Analysis 

Figure A-1 is the block diagram of the closed-loop, noise-injection, 
feedback deterministic radiometer model upon which this analysis is 
based (ref. 1). For analysis purposes, this block diagram was simplified 
to the block diagram of figure A-2 in which Boltzmann’s constant, the 
square-law detector constant, the gains of the RF and AF amplifiers, the 
Dicke constant, and the gains of the integrator and the lead- lag network 
have all been combined into Ki, and the gains of the directional 
coupler and isolator and the conversion factors associated with the noise 
diode have been combined into K 2 . 

A basic measure of system performance used in the evaluation was 
the response of the system to a unit step input. If such an input is 
applied to the system of figure A-2, some steady-state output will 
eventually be reached as indicated in figure A-3. If we define settling 
time, t^, as the time required for the output to reach and remain 
within some arbitrary percentage of the final value, settling time can be 
used as a measure of system performance. For this evaluation, the 
settling time was defined as the time after which the output would not 
deviate from the final value by more than 0.01 percent. 

The noise-equivalent bandwidth, of a system with transfer 

function H(f), is the bandwidth of an ideal filter with the same midband 
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gain, as H(£), and which passes the same average noise power 

as HCf) when white noise is applied as its input. The standard 
definition of B. is 


B 


N 


o 


/ 


oo 


|HCf)|2 


df 


(A-1) 


To determine the noise-equivalent bandwidth of the system of figure A-2, 
temporarily use p as the independent variable of the transfer function. 
Then, H(p) is given by 


H(p;) 


1 + 


(1 + T P) 

K2 2 

mcTKifT^ tT" 
K^Ki /P * 


P 


2 


1 

K2 


+ 0) T 
o 


2 



1 




(A-2) 


where 



1 + KJK2T2 


and 

• o 


K1K2 


Making the change of variable, s = 


H(s) = (A-3) 

1 + 2cs + s^ 
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Evaluating equation (A-1) for this transfer function 


B 


N 


1 + (o)^T2)2 

' 8c 


(A-4) 


where is a measure of the average noise power that will be present 

in the output of the system when white noise is present at the input. 

In general, as the parameters of the system of figure A-2 are 
varied, both B^^ and t^ will change. Although optimum performance 
involves minimizing both B^^ and t^, parameter values which minimize 
Bj^ tend to make t^ large, and vice versa. Therefore the product of 
Bj^ and t^ was selected as a performance criterion. 

The underdamped, critically damped, and overdamped responses to 
the system of figure A-2 to a unit step input are given by equations 
(A-5), (A-6) and (A-7) respectively. In these equations, the actual 
unit step response is yi (t) where y(t) = K 2 yi (t) . 

Underdamped response 


y(t) = 1 


Ae 


-?t 


sin 




t + 0 


(A-5) 


where 


A 



+ 1 


and 


0 


tan 


Vi - ) 

c - 0)^T2/ 
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I Mil 


IIIH Hill 


Critically damped response 


y(t) = 1 - [(1 - (o^T 2 )t + l] e 
Overdamped response 


y(t) 


1 


Be sinh 



1 t + (|) 


(A-6) 


(A- 73 


where 


B 



C? - 

” 7 ^' -""i “ 


2 


- 1 


and 


- “o^2 

Figure A-4 indicates how the settling time was determined for the 
underdamped case where the overshoot exceeds the target percentage of 
deviation (figure A-4a3, the underdamped case where the overshoot is less 
than the target percentage of deviation (figure A-4b) , and the critically 
damped and overdamped cases (figure A-4c). 

Since both B^ and t^ are functions of the damping factor, c, 
changes in the damping factor can be expected to change the 
product. If the per unit sensitivity of damping factor to is 

defined as 


A? 




(A-83 
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where Kj^ = K 1 K 2 , the sensitivity of damping factor to changes in loop 
gain for the system of figure A- 2 is 


S 


1 1 

2 1 + Kj^T2 


CA-9) 


which can be expressed as 


S 





u) To 
o 


(A- 10) 


Maximum overshoot and the time at which it occurs for the under- 
damped case were evaluated by setting the partial derivative of y(t) 
with respect to t equal to zero and solving the resulting equation 
for t. The result is equation (A-11). This value of t is substituted 
into the expression for l-y(t) to obtain the maximum overshoot of 
equation (A- 17). 


max 


Vi - ?■ 


tan 


-1 


'^2- 1 


+ nil 


(A-11) 


0 Cw To > 1 
o — 


1 ?U) To < 1 

o ^ — 


max 


= I cos ^ 


1 - t + 


max 


Vl - 


Sin 


t \ 
y max [ 


-?t 


max 


(A-12) 


Defining the error as l-y(t), the per unit sensitivity of the 
error to changes in Kj can be evaluated. Define the per unit sensitivity 
of the error to changes in as 
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(A- 13) 
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After making the substitution s = jw, this transfer function can be 
put in the following form: 


1 + j 2 — _ 

1 - 

The half-power points for this function occur when 

2^(f/f ) 

= ±1 (A-19) 

1 - 

so that for positive frequencies, the half-power points are at the 
following frequencies: 

£2 = £„(c * i) 

£l = £j,(-C + l) (A-20) 

the bandwidth of the filter, BW = £2 - fi, is then 

BW = 2cf^ (A-21) 

The response of this filter to a unit step input is 


y(t) = 1 




e 


t 

o 


for the underdamped case. 


(A-22) 
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Let 


e 



sin 




-Ca.ot 


CA-23) 


be the error. For the underdamped case this error will have a 
maximum value when 


t = tan ^ 



CA-24) 


as indicated in figure A-5a, after which the error will have a series 
of peaks of alternating sign and decreasing magnitude. Consider the 
first negative maximum of the error. Because of the exponential factor 
in the expression for the error, no succeeding maxima or minima will 
exceed it in magnitude, so it represents the maximum overshoot of the 
step response, the maximum amount by which the output exceeds the final 
value of the output, although it is not the largest error. In figure 
A-5a, the settling time tj is the time after which the error never 
exceeds the maximum overshoot. 

For the critically damped and overdamped cases, the response is 
as indicated in figure A-5b. 

PRESENTATION OF RESULTS 


B.,t Product 
N s 

In figure A-6, the product of the noise-equivalent bandwidth B^^ 
and the settling time t^ is plotted versus the damping factor ? for 
five values of lead time constant. The simple lag network, in which 
the lead time constant T 2 is zero, appears best from a consideration of 
just the product, since it not only has the lowest values, but 
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the values remain relatively constant over a wide range of damping factor. 

The abrupt decreases in the product that are apparent in the 

curves for o) To = 2.0 and co To = 1.5 occur when 
o o ^ 


C = 


1 + Chi^T2) 
2io^T2 


Sensitivity of Damping Factor to Kj^ Changes 

The sensitivity of damping factor to changes in the loop gain is 
plotted as a function of damping factor for five values of lead time 
constant in figure A-7. 


Sensitivity of Error to Changes 

Figure A- 8 is a plot of the sensitivity of the error to changes in 
the forward gain plotted as a function of time for three lead time 
constants and for a damping factor of 0.95, the value at which the noise- 
equivalent bandwidth and settling time product is minimum when the lead 
time constant is zero. 


Damping Factor and Settling Time versus Maximum Overshoot 

Figures A-9, A-10, A-11, A-12, and A-13 are plots of damping factor 
and versus maximum overshoot for the underdamped loop, each 

figure for a different lead time constant. A comparison of the effect 
of lead time constant on overshoot and the resultant settling time can 
be obtained from these curves. 

Notch Filter Damping Factor and Settling Time 

Figure A- 14 presents damping factor and settling time as a function 
of maximum overshoot for the notch filter. 
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SUMMARY 


Loop Analysis 

Figures A-6, A-7, and A-8 indicate that as long as loop gain vari- 
ations are not encountered, the simple lag network provides superior 
performance as long as damping factors above approximately 0.8 are 
used. For damping factors in the vicinity of critical damping, which 
appears to be the logical area of damping factor to operate, = 1.0 

minimizes the sensitivity of the damping factor to changes in loop gain. 
Although figure A-8 indicates that the error is still very sensitive to 
changes in Kj at most times, this is probably not as serious as it 
appears at first. 

For example, considering the curve of figure A-8 for = 1.0, 

it can be seen that the error is insensitive to changes in Kj in the 
vicinity of w^t = 10.5. From figure A- 11 it can be seen that the 
maximum overshoot corresponding to a settling time of this value is 
less than 10"^. Since can be expected to be on the order of 100 

rad/s, this operation would require allowing times on the order of 
0.1 s for settling, which is not unreasonable. 

Although figure A-6, which is plotted for an error of 10 indicates 
that the simple lag network produces the minimum bandwidth-time product 
for any damping factor, the sensitivity of the error to changes of 
loop gain will likely preclude operating the system at minimum settling 
times. Allowing additional time for settling before taking a measure- 
ment will reduce the magnitude of the error. However, so much time may 
be required that lead compensation may be preferable, even though its 
minimum bandwidth-time product is higher than that of the simple lag 
network. 
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Figure A-1. Block diagram of closed-loop, noise-injection, feedback deterministic 
radiometer models. 
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Figure A- 14. Damping factor and settling time versus maximum 
overshoot for the notch filter. 
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TABLE 1.- PREDICTED SENSITIVITIES BASED ON NOISE FLUCTUATIONS AND MEASUREMENT TIMES 
(INCLUDING BOTH PROCESSING AND SETTLING TIMES) FOR PROPOSED DIGITAL SYSTEM. 





N=1 


2 


4 


8 


16 


32 


64 




T 

26.75 

ms 

53.5 

ms 

107 

ms 

214 

ms 

428 

ms 

856 

ms 

1.712 

s 



AT 

2.096 

K 

1.504 

K 

1.168 

K 

0.866 

K 

0.626 

K 

0.447 

K 

0.318 

K 


B 
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■"s 

146 

ms 

153 

ms 

158 

ms 

164 

ms 

170 

. ms 

177 

ms 

184 

ms 

= 20 

MHz 



















173 

ms 

206.5 

ms 

265 

ms 

378 

ms 

598 

ms 

1.033 

s 

1.896 

s 



AT 

0.958 

K 

0.672 

K 

0.523 

K 

0.387 

K 

0.280 

K 

0.200 

K 

0.142 

K 


B 

















SI 
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162 

ms 

168 

ms 

174 

ms 

180 

ms 

186 

ms 

193 

ms 

200 

ms 

= 100 

MHz 


















"m 

189 

ms 

221.5 

ms 

281 

ms 

394 

ms 

614 

ms 

1.049 

s 

1.912 

s 



AT 

0.419 

K 

0.301 

K 

0.234 

K 

0.173 

K 

0.125 

K 

0.0894 

K 

0.0636 

K 


B 

















SI 


















■^s 

178 

ms 

185 

ms 

190 

ms 

196 

ms 

202 

ms 

209 

ms 

216 

ms 

= 500 

MHz ! 


















"m 

205 

ms 

238.5 

ms 

297 

ms 

410 

ms 

630 

ms 

1.065 

s 

1.928 

s 



AT 

0.210 

K 

0.150 

K 

0.117 

K 

0.0866 

K 

0.0626 

K 

0.0447 

K 

0.0318 

K 
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■^s 

192 

ms 

199 

ms 

204 

ms 

210 

ms 

216 

ms 

223 

ms 

230 

ms 

= 2 

GHz 



















219 

ms 

252.5 

ms 

311 

ms 

424 

ms 

644 

ms 

1.079 

s 

1.942 
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TABLE 2.- COMPUTER SIMULATION OF PROPOSED DESIGN WITH TWO POST-LOOP SAMPLES 




IN AVERAGE. (ALL UNITS 

EXCEPT TIME ARE EXPRESSED IN 

KELVIN.) 


72 SIMULATION 

TIME(s) 



Temperature Estimates 


y 

2 

s 

^2 

a 

^2 

Run 1 

Run 2 

Run 3 

Run 4 

0 

213 

216 

54 

188 

.11 

188 

29 

248 

.71 

210 

.41 

28 

.82 

24 

.06 

0 

.320 

193 

82 

180 

90 

243 

18 

190 

.79 

202 

.17 

27 

.89 

24 

.06 

0 

.427 

203. 

17 

221 

.00 

196 

04 

172 

.97 

198 

.30 

19 

.88 

24 

.06 

0 

.533 

166 

72 

218 

59 

213 

69 

177 

51 

194 

13 

25 

.87 

24 

.06 

0 

.640 

162. 

72 

201 

21 

210 

12 

200 

32 

193 

.59 

21 

.05 

24 

.06 

0 

747 

212. 

08 

177 

42 

207 

01 

185 

71 

195 

56 

16 

63 

24 

06 

0 

853 

196. 

40 

225 

45 

236. 

77 

147 

30 

201 

48 

39 

92 

24 

.06 

0 

960 

257. 

53 

222 

15 

212. 

35 

223 

05 

228 

77 

19 

78 

24 

06 

1 

067 

245, 

77 

222 

87 

186, 

42 

229 

46 

221 

13 

25 

06 

24 

06 

1 

173 

158. 

80 

197. 

65 

238. 

19 

190 

61 

196. 

31 

32 

64 

24. 

06 

1. 

280 

202. 

10 

213. 

07 

188. 

74 

201. 

84 

201, 

44 

9 

95 

24 

06 

I. 

387 

202. 

28 

195. 

78 

219. 

39 

224 

02 

210. 

37 

13. 

49 

24. 

06 

1. 

493 

196. 

76 

269. 

65 

216. 

54 

197. 

65 

220. 

15 

34. 

24 

24. 

06 

1. 

600 

212. 

08 

178. 

22 

216. 

90 

213. 

33 

205, 

13 

18. 

06 

24. 

06 

1. 

707 

212. 

80 

169. 

85 

170. 

65 

209. 

95 

190. 

81 

23. 

77 

24. 

06 

1. 

813 

198. 

18 

212. 

00 

180. 

45 

161. 

74 

188. 

09 

21. 

80 

24. 

06 

1. 

920 

184. 

01 

239. 

98 

177. 

06 

214. 

31 

203. 

84 

29. 

02 

24. 

06 











203. 

63 

23. 
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TABLE 3.- COMPUTER SIMULATION OF PROPOSED DESIGN WITH FOUR POST-LOOP SAMPLES 


IN AVERAGE. (ALL UNITS EXCEPT TIME ARE EXPRESSED IN KELVIN.) 


SIMULATION 



Temperature 

Estimates 





TIME(s) 

Run 1 

Run 2 

Run 3 

Run 4 


s 

^74 

0.213 

202.59 

211.73 

208.70 

228.44 

212.87 

11.06 

18.69 

0.320 

196.85 

179.48 

240.51 

171.67 

197.13 

30.78 

18.69 

0.427 

198.09 

241.05 

194.89 

184.59 

204.66 

24.94 

18.69 

0.533 

171.00 

222.24 

206.20 

174.57 

194.00 

25.62 

18.69 

0.640 

184.59 

195.02 

188.20 

206.83 

193.50 

24.85 

18.69 

0.747 

211.10 

186. 15 

205.45 

180.54 

195.81 

14.76 

18.69 

0.853 

200.37 

205.67 

217.48 

169.31 

198.21 

20.55 

18.69 

0.960 

229.19 

216.50 

205.67 

204.15 

213.88 

11.60 

18.69 

1.067 

231.69 

224.20 

192.26 

224.74 

218.22 

17.64 

18.69 

1.173 

153.94 

211.19 

234.36 

194.80 

198.57 

33.89 

18.69 

1.280 

199.70 

208.21 

208.34 

186.33 

200.65 

10.36 

18.69 

1.387 

203.80 

199.74 

216.81 

201.93 

205.57 

7.67 

18.69 

1.493 

214.71 

217.48 

204.91 

207.01 

211.03 

6.02 

18.69 

1.600 

191.23 

217.83 

217.43 

217.83 

211.08 

13.23 

18.69 

1.707 

190.83 

174.44 

188.07 

211.64 

191.25 

15.37 

18.69 

1.813 

195.24 

211.42 

203.93 

192.84 

200.86 

8.50 

18.69 

1.920 

188.96 

235.52 

178.62 

207.90 

202.75 

24.99 

18.69 






202.94 

17.75 







A 

^4 

A 

Sy 

•^4 



£ 




TABLE 4.- COMPUTER SIMULATION OF PROPOSED DESIGN WITH EIGHT POST-LOOP SAMPLES 
IN AVERAGE. (ALL UNITS EXCEPT TIME ARE EXPRESSED IN KELVIN.) 


78 SIMULATION 



Temperature 

Estimates 





TIME(s) 

Run 1 

Run 2 

Run 3 

Run 4 

^8 

^>^8 . 

0 

^8 

0.320 

199.72 

195.60 

224.61 

200.05 

205.00 

13.23 

13.86 

0.427 

197.47 

210.26 

217.70 

178.13 

200.89 

17.32 

13.86 

0.533 

184.55 

231.64 

200.54 

179.58 

199.08 

23.48 

13.86 

0.640 

177.80 

208.63 

197.20 

190.70 

193.58 

12.87 

13.86 

0.747 

197.85 

190.59 

196.82 

193.68 

194.74 

3.28 

13.86 

0.853 

205.74 

195.91 

211.46 

174.93 

197.01 

16.06 

13.86 

0.960 

214.78 

211.08 

211.57 

186.73 

206.04 

12.98 

13.86 

1.067 

230.44 

220. 35 

198.96 

214.45 

216.05 

13.17 

13.86 

1.173 

192.81 

217.70 

213.31 

209.77 

208.40 

10.89 

13.86 

1.280 

176.82 

209.70 

221.35 

190.56 

199.64 

19.80 

13.86 

1.387 

201.75 

203.98 

212.58 

194.13 

203.11 

7.59 

13.86 

1.493 

209.26 

208.61 

210.86 

204.47 

208.30 

2.72 

13.86 

1.600 

202.97 

217.65 

211.17 

212.42 

211.05 

6.08 

13.86 

1.707 

191.03 

196. 13 

202.75 

214.74 

201.16 

10.24 

13.86 

1.813 

193.04 

192.93 

196.00 

202.24 

196.05 

4.36 

13.86 

1.920 

192.10 

223.47 

191.28 

200.37 

201.81 

15.02 

13.86 






202.62 

11.86 


1 





A 

^8 

A 

s 

>^8 
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TABLE 5.- COMPUTER SIMULATION OF PROPOSED DESIGN WITH 16 POST-LOOP SAMPLES IN 
AVERAGE. (ALL UNITS EXCEPT TIME ARE EXPRESSED IN KELVIN.) 


yi6 SIMULATION 



Temperature 

Estimates 





TIME(s) 

Run 1 

Run 2 

Run 3 

Run 4 


s 

^16 

a 

^16 

0.533 

192.14 

213.62 

212.58 

189.82 

202.04 

12,81 

10.02 

0.640 

187,64 

209.44 

207.45 

184.42 

197.24 

13.03 

10.02 

0.747 

191.20 

211.12 

198.68 

186.63 

196.91 

10.70 

10.02 

0. 853 

191.77 

202.27 

204.33 

182.81 

195.30 

9.98 

10.02 

0.960 

206.32 

200.83 

204.20 

190.21 

200.39 

7.15 

10.02 

1.067 

218.09 

208.13 

205.21 

194.69 

206.53 

9.63 

10.02 

1.173 

203.80 

214.39 

212.44 

198.25 

207.22 

7.55 

10.02 

1.280 

203.63 

215.03 

210.16 

202.51 

207.83 

5.87 

10.02 

1.387 

197.28 

210.84 

212,94 

201.95 

205.75 

7.39 

10.02 

1.493 

193.04 

209.16 

216.11 

197.52 

203.96 

10.57 

10.02 

1.600 

202.36 

210.82 

211.87 

203.27 

207.08 

4.96 

10.02 

1.707 

200.14 

202.37 

206.81 

209.60 

204.73 

4.27 

10.02 

1.813 

198.01 

205.29 

203.59 

207.33 

203.56 

4.00 

10.02 

1.920 

191.57 

209.80 

197.01 

207.55 

201.48 

8.65 

10.02 






202,86 

A 

^16 

8.33 

^16 





TABLE 6.- COMPUTER SIMULATION OF PROPOSED DESIGN WITH 32 POST-LOOP SAMPLES 
IN AVERAGE. (ALL UNITS EXCEPT TIME ARE EXPRESSED IN KELVIN.) 


732 SIMULATION 



Temperature Estimates 





TIME(s) 

Run 1 

Run 2 

Run 3 

Run 4 

y 32 

^32 

a 

^32 

0.960 

199.23 

207.23 

208.39 

190.01 

201.22 

8.51 

7.15 

1.067 

202.86 

208.79 

206.33 

189.55 

201.88 

8.57 

7. 15 

1.173 

197.50 

212.75 

205.56 

192.44 

202.06 

8.94 

7.15 

1.280 

197.70 

208.65 

207.25 

192.66 

201.57 

7.68 

7.15 

1.387 

201.80 

205.83 

208.57 

196.08 

203.07 

5.43 

7.15 

1.493 

205.56 

208.64 

210.66 

196.10 

205.24 

6.44 

7.15 

1.600 

203.68 

212.60 

212.16 

200.76 

207.30 

5.99 

7.15 

1.707 

201.89 

208.70 

208.48 

206.05 

206.28 

3.16 

7.15 

1.813 

197.64 

208.06 

208.26 

204.64 

204.65 

4.96 

7.15 

1.920 

192.30 

209.48 

206.56 

202.53 

202.72 

7.51 

7.15 






203.60 

6.72 







A 

>^32 

A 

Sy 
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TABLE 7.- COMPUTER SIMULATION OF PROPOSED DESIGN WITH 64 POST-LOOP SAMPLES IN 
AVERAGE. (ALL UNITS EXCEPT TIME ARE EXPRESSED IN KELVIN.) 




y 

04 SIMULATION 






Temperature 

Estimates 





TIME(s) 

Run 1 

Run 2 

Run 3 

Run 4 

^64 

^^64 

'"y64 

1.813 

198.43 

207.65 

208.33 

196.04 

202.61 

6.29 

5.09 

1.920 

197.58 

209.13 

206.45 

197.32 

202.62 

6.07 

5.09 






202.62 

6.18 


1 

1 





A 

^64 

A 

•^64 
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1 = voltage- to- frequency 

^ ^ / pulses/s\ 

constant I *- — ^ j 

= pulse width (s/pulse) 



1 


voltage-to- frequency constant 


pulses/s 

V 



DUTY 

^CYCLE 


nx 


T 


= pulse width (s/pulse) 

*DC>V ♦ * ■'n) 


E ■ !<’' - V- \ - 
N = WeX^'^'"o'= 

^ = Sc'^ATT^EX^'^l 


Vo (l-Kj^)G/2 

^B'^A l+GB/2 

Fo (l-Kp^)GKi/2 

Tg-T^ 1 + G 6/2 

•^E 

T^-Tp 1 + G6/2 


1 


K 

1 




Figure 1. Block diagram of control mechanism in feedback noise-injection radiometer 











107 


1 GHz 



I 

I 


Figure 2. Block diagram of radiometer system used as a reference in this study. 













30 KU 


C2 



Rj = 

Cl = 1 mF 
Rz~ 500 Kfi 
C 2 = 0.001 viF 


b) AMPLITUDE RESPONSE 



25.95 dB (19.8') 


5628 Hz 



Figure 4. Circuit used for stage IC3 and its Bode plot. 
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0-300 K 



Figure 9. Illustration of noise addition in loop and range 
of values. 
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Figure 11. Control loop model for stepped frequency radiometer 
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Figure 19. Sampling of bandwidth- limited signal exactly at Nyquist 

rate (f = 2 B. ) . 
s 1 
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Figure 20. Sampling of bandwidth- limited signal at one-fourth the 
Nyquist rate. 


126 




B. 

1 



f -B. 
s 1 


Figure 21. Sampling of bandwidth- limited signal at twice the Nyquist 
rate. 


127 







Figure 22. Sampling at different rates of signal whose spectrum has 
gradual rolloff. 
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Figure 24. Mode A of the second second-order digital estimating circuit investigated. 
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Figure 28, Proposed loop estimation and feedback system 

















T = 365,000 K 


Figure 30. Ranges of temperatures and duty cycle of proposed design. 
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Figure 32. Block diagram of proposed system 


HOT CARRIER DIODE 



J 

POST- LOOP 

n 

n 

PROCESSOR 


DATA 

OUT 




CONTROI. 

FEEDBACK 
















139 



■riJE ACTUAL LOOP OUTPUT y(n) CHANGES EVERY 1.6718 ms. 
HOWEVER, ALL yjj(i) ESTIMATES CHANGE EVERY 26.75 ms. 



NUMBER 

ESTIMATE OF POINTS 




ye(i) 

XieCi) 

Y32Ci) 

y6**(i) 



I 


2 


4 


8 


16 


32 


64 


for post- loop processor. 
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Figure 34. Digital filter format for sum- and -dump algorithms. 





















Figure 36. Sensitivities of proposed design resulting from dynamic 
simulation (78.125-kHz bandwidth). 
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